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A.  Scientific  Research  Goals 

Precipitation  strengthening  is  a  potentially  attractive  technique  for 
hardening  and  toughening  infrared  transmitting  materials.  The 
controlled  precipitation  of  a  finely  dispersed  second  phase  offers  the 
opportunity  to  improve  mechanical  properties  without  degrading  optical 
properties.  It  is  well  established  that  the  presence  of  an  appropriate 
dispersion  of  second  phase  particles  can  improve  the  mechanical  behavior 
of  ceramics,  and,  by  maintaining  a  particle  size  <100  nm,  scattering  losses 
in  the  infrared  will  be  low  and  the  need  to  index-match  the  ZnS  matrix  with 
the  second  phase  will  be  circumvented. 

The  overall  objectives  of  the  research  program  have  been  (1)  to 
demonstrate  that  ZnG-based  ceramics  can  be  precipitation-strengthened 
and  (2)  to  characterize  the  resulting  mechanical  and  optical  properties  of 
the  multiphase  materials.  An  important  research  goal  has  been  to  identify 
ZnS-based  systems  suitable  for  strengthening  and  to  determine  the 
composition  ranges  and  annealing  treatments  required  to  produce 
precipitate  phases  or  other  solid-state  reaction  products.  Another 
important  research  objective  has  been  to  determine  why  "alloys"  in  the  ZnS- 
Ga2S3  system  exhibit  an  increase  in  both  fracture  toughness  and  hardness, 
and  to  establish  the  composition  ranges  and  annealing  treatments  which 
produce  this  behavior.  The  research  has  also  sought  to  quantify  the  effects 
of  size  and  volume  fraction  of  various  second  phase  constituents  on  infrared 
scattering.  The  relationship  between  mechanical  behavior  and 
microstructural  development  has  been  of  central  importance  in  this 
research  program. 

B.  Summary  of  Research  Accomplishments 

Candidate  systems  for  precipitation  strengthening  have  been 
identified,  i.e.,  the  ZnS-CdS  and  ZnS-ZnGa2S4  systems.  The  phase 
equilibria  in  the  ZnS-rich  region  of  these  two  systems  have  been  established. 
CdS  is  very  soluble  in  ZnS  and  the  two-phase  microstructure  envisioned  in 
this  system  will  consist  of  sphalerite  and  wurtzite  solid  solutions.  Ga2S3 
exhibits  limited  solubility  in  ZnS  and  the  pseudobinary  ZnS-ZnGa2S4  system 
is  characterized  by  a  eutectoid  reaction.  The  equilibrium  phase  diagrams 
for  both  systems  were  published  and  ar^  attached  to  this  report. 


The  research  program  focussed  cn  the  ZnS-ZnGa2S4  system.  We 
developed  hot  pressing  methods  for  ZnS  and  ZnS-Ga2S3  "alloys"  which  4 
routinely  produced  samples  exceeding  99.0%  of  their  theoretical  density. 
These  high  quality  specimens  were  used  to  characterize  mechanical  and 
optical  properties  of  materials  containing  8  to  16  mol%  Ga2S3. 

The  most  significant  result  obtained  during  the  program  was  the 
observation  that  the  various  ZnS-Ga2S3  compositions  exhibit  considerably 
higher  fracture  toughness  and  hardness  than  pure  ZnS.  The  increase  in 
hardness  was  proportional  to  the  amount  of  Ga2S3  added,  with  the  16  mol% 
sample  exhibiting  a  value  nearly  two  times  larger  than  ZnS.  A  sizeable 
increase  in  fracture  toughness  was  also  observed  in  all  samples.  At  12 
mol%  Ga2S3,  the  fracture  toughness  value  was  more  than  50%  greater  than 
that  of  pure  ZnS.  In  these  studies,  all  samples  were  solid-solution  treated 
in  the  wurtzite  phase  (at  940°C),  and  no  precipitated  phases  were  detected  by 
x-ray  diffraction.  The  mechanism  responsible  for  the  observed  increase  in 
hardness  and  fracture  toughness  has  yet  to  be  determined.  The  publication 
describing  these  results  is  attached  to  the  report. 

The  optical  properties  of  several  ZnS-Ga2S3  solid  solutions  have  been 
measured.  As  expected,  the  infrared  transmission  exhibited  no  change 
from  that  observed  for  pure  ZnS,  with  excellent  transmission  until  the 
infrared  cut-off  in  the  12  pm  range.  In  addition  to  this  experimental  work, 
a  model  has  been  derived  which  enables  us  to  calculate  the  scattering 
efficiency  factor  as  a  function  of  wavelength  for  a  multiphase  material  with 
a  given  particle  size  and  refractive  index  for  the  second  phase.  By  applying 
this  model  to  pure  ZnS  (where  pores  constitute  the  second  phase),  we  were 
able  to  make  a  unique  calculation  of  the  pore  size  distribution  in  the 
ceramic  by  measuring  the  optical  transmission.  A  manuscript  describing 
the  model  is  attached  to  this  report. 

Another  significant  result  is  the  successful  development  of  a  method 
far  measuring  fracture  toughness  using  an  indented  miniaturized  ring-on- 
ring  disk-bend  fracture  toughness  test  (MDBFTT).  Samples  3  mm  in 
diameter  and  300  to  400  pm  thick  are  indented  using  a  Vickers  indenter 
and  the  load  to  fracture  is  measured  in  the  MDBFTT.  This  is  repeated  for 
several  loads,  af^r  whirh  ~lrp ady  Hovplnnpd  fr^rturp  toughness  equations 
are  applied  to  the  data,  yielding  a  value  of  Kic.  The  fracture  toughness  can 
be  obtained  by  this  method  without  measuring  actual  crack  lengths,  thereby 


eliminating  the  uncertainties  associated  with  such  measurements.  Initial 
measurements  on  standard  samples,  including  CVD-grown  ZnS,  have 
yielded  values  which  are  in  excellent  agreement  with  conventional  fracture 
toughness  measurements.  The  attached  manuscript  provides  full  details  of 
the  technique  and  the  results  to  date. 
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The  ZnS-GajS.i  equilibrium  phase  diagram  has  been  deter¬ 
mined  to  511  mol'?}-  over  the  temperature  range  700°  to 
900°C.  Samples  of  various  compositions  were  prepared  via 
solid-state  diffusion  starting  from  powders  of  the  pure  com¬ 
ponents.  The  identification  of  the  phases  was  determined  by 
X-ray  diffraction  methods.  The  principal  feature  of  the 
phase  equilibria  is  the  eutectoid  transformation  at  818  ±  5°C 
of  hexagonal  wurtzite  containing  16  ±  1  mol%  Ga2Sj  to 
..uhic  7.nS  and  tetragonal  ZnGa2S4.  ZnGa2S4  is  the  equi¬ 
librium  compound  at  50  mol%  Ga.S.(,  but  it  exists  over  a 
considerable  range  of  stoichiometry.  The  solubility  of  Ga2S2 
in  ZnS  increases  with  increasing  temperature  to  a  maxi- 
mumofO  ±  1  mol%  at  the  eutectoid  temperature.  [Key  words: 
phase  diagrams,  zinc  sulfide,  solid  state,  wurtzite,  eutectics.| 

I.  Introduction 

Zinc  si  i  i  idi  is  commercially  used  as  an  infrared  (IR) 
transmitting  material  because  of  its  good  optical  trans¬ 
mission  up  to  12  p im  and  adequate  mechanical  strength.  How¬ 
ever.  ZnS-hased  ceramics  with  improved  fracture  toughness, 
thermal-shock  resistance,  and  erosion  resistance  are  abso¬ 
lutely  essential  for  certain  IR-window  applications.  There  are 
several  possible  approaches  to  attain  the  desired  combination 
of  improved  mechanical  properties  while  retaining  IR  trans¬ 
mittance  comparable  to  that  of  ZnS.  Our  approach  is  to  ex¬ 
plore  the  possibilities  of  strengthening  and  toughening  ZnS  by 
the  addition  of  dispersions  of  particles  of  a  second  phase,  with 
the  dispersoids  introduced  into  the  microstructure  by  the  ap¬ 
propriate  heat  'reatment  of  densified  samples.  The  second- 
phase  additions  should  also  have  intrinsically  good  IR 
transmittance  and  a  closely  matching  refractive  index  in  order 
to  minimize  optical  attenuation  losses  due  to  scattering.  This 
requirement  places  certain  restrictions  on  the  suitability  of  the 
ions  added  for  alloying  purposes. 

To  improve  the  mechanical  behavior  of  ZnS  through  heat 
treatment,  knowledge  of  solid-solid  phase  equilibria  involving 
ZnS  and  other  sulfides  is  also  essential.  Unfortunately,  such 
data  are  limited,  and  we  have  had  to  undertake  the  task  of 
determining  some  binary  phase  diagrams  ourselves.  Among 
the  various  possibilities  we  have  considered,  the  ZnS-C’dS 
and  ZnS-Ga.>Si  systems  appear  to  be  especially  promising  for 
meeting  all  the  requirements  noted  above.  We  have  already 
reported  the  results  of  the  determination  of  the  ZnS-rich  por¬ 
tion  ol  the  ZnS-CdS  phase  diagram1  and  have  presented  a 
preliminary  version  of  the  ZnS-Cia.-Si  phase  diagram.'  In  this 
paper,  we  present  the  ZnS-Cia:Si  phase  diagram  to  50  moUf 
Ga  .  N  i.  at  'he  same  time  refining  the  earlier  version. 

I  he  first  study  of  this  phase  diagram  was  made  by  Gates 
and  Edwards,'  using  the  simultaneous  Knudson  anti  dynamic 
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torsion-effusion  method.  They  published  the  Ga-Si-rieh  end 
of  the  partial  ZnS-Ga2Sv  phase  diagram,  locating  the  phase 
boundaries  to  within  ±15  C.  However,  this  part  of  the  dia¬ 
gram  is  of  no  interest  for  strengthening  ZnS-based  ceramics. 
Malevskii4  studied  the  ZnS-Ga;S4  system  using  the  solid- 
state  reaction  technique  and  published  a  partial  phase  dia¬ 
gram  at  the  ZnS-rich  ,_nd.  His  work  has  established  that  the 
solubility  of  Ga-Si  increases  with  increasing  temperature  and 
that  the  diagram  contains  a  eutectoid  reaction  around  SIX)  C. 
The  eutectoid  reaction  reported  involved  the  omposition 
of  wurtzite  containing  approximately  20  molrz  Ga2S-.  to  form 
sphalerite  and  the  thiogailate  phase  ZnGa:S4.  containing  ap¬ 
proximately  14  and  2<S  molG  Ga-Si.  respectively. 

ZnGa;Sj  itself  is  a  good  IR-transmitting  material  in  the 
range  2.5  to  12  gm  and  has  the  742 m  tetragonal  crystal  struc¬ 
ture.'''  The  tetragonal  unit  cell  has  lattice  parameters  a,  and 
c,  of  0.5297  and  1.0363  nm,  respectively.1  The  lattice  mismatch 
with  cubic  sphalerite  ZnS  (a,  =  0  5410  nm)  is  very  small, 
which  probably  accounts  for  the  fairly  large  solubility  of 
thiogailate  in  sphalerite  (despite  the  difference  in  valence  be¬ 
tween  Zn  and  Ga)  reported  by  Malevskii.  Also,  the  micro- 
hardness  of  ZnGa2S4  has  been  found  to  be  close  in  value  to 
that  of  ZnS.7 

Based  on  the  factors  discussed  above,  it  is  clear  that  Ga;S> 
satisfies  the  requirements  for  a  suitable  second-phase  addi¬ 
tion.  ZnS-rich  solid  solutions  should  be  precipitation  harden- 
able.  whereas  alloys  of  near  eutectoid  composition  offer  the 
prospect  of  microstructural  manipulation  by  eutectoid  de¬ 
composition.  Additionally,  the  expected  decomposition 
product.  ZnGa.-Sj,  is  optically  compatible  with  ZnS.  One  of 
the  major  drawbacks  to  the  development  of  these  "alloys"  as 
IR  materials  is  the  accuracy  of  Malevskiis  phase  diagram, 
because  the  compositions  of  the  samples  used  in  his  investiga¬ 
tion  differed  by  as  much  as  10  molT.  This  was  the  main  rea¬ 
son  for  undertaking  the  current  investigation. 

II.  Experimental  Procedure 

X-ray  diffraction  was  the  principal  technique  employed  to 
investigate  the  phase  equilibria.  Powders  of  ZnS  (40.00  )*  and 
Ga.-Si  (49. WO)f  were  used  as  starting  materials  for  the  prepa¬ 
ration  ol  alloys  by  solid-state  reaction.  The  method  of  prepa¬ 
ration  used  in  the  present  study  was  previously  described  by 
Chen  ct  al.'  The  powders  were  mixed,  ground  with  a  mortar 
and  pestle,  sealed  under  vacuum  in  a  fused  silica  ampule,  and 
heated  to  the  reaction  temperature.  Independent  chemical 
analyses  were  not  performed  on  the  reacted  powders,  so  the 
concentrations  reported  are  those  derived  from  the  original 
measurements  of  the  weighed  powders.  The  accuracy  of  these 
measurements  is  ±0.5'~; . 

In  order  to  assure  that  equilibrium  was  attained,  the  reac¬ 
tion  times  were  chosen  so  that  the  lattice  constants  of  the 
solid  solutions  no  longer  changed  as  the  reaction  lime  in¬ 
creased.  A  solid-state  reaction  period  of  more  than  It)  weeks 
was  necessary  at  700  (',  but  shorter  times  were  used  at  the 
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higher  reaction  temperatures.  The  ampules  were  quenched 
into  water  at  room  temperature  after  all  reaction  treatments. 
The  various  conditions  are  summarized  in  Table  I.  The  inter¬ 
vals  in  composition  of  the  samples  were  1  or  2  mol%  near  the 
pha'e  boundaries  and  3  or  5  mol%  elsewhere.  We  have  found 
that  the  approach  to  equilibrium  is  far  fastu  for  the  mixed 
powders  than  for  compacted  samples  of  identical  composi¬ 
tion.  This  is  because  the  eutectoid  decomposition  reaction  is 
very  sluggish. 

The  parametric  method"  and  a  modified  version  of  the 
disappearing-phase  method"1'  were  used  to  establish  the  solvus 
compositions.  The  choice  of  method  was  determined  by  the 
crystal  structures  of  the  equilibrium  phases  present.  In  par¬ 
ticular,  the  X-ray  diffraction  peaks  of  wurtzite  are  readily 
distinguished  from  those  of  the  sphalerite  and  thiogaliate 
phases.  However,  the  sphalerite  and  thiogaliate  phases  are 
structurally  similar,  and  many  of  their  peaks  overlap.  For  this 
reason,  the  disappearing-phase  method  was  used  to  deter¬ 
mine  the  phase  boundaries  at  850°  and  900°C,  where  wurtzite 
is  in  equilibrium  with  either  sphalerite  or  thiogaliate,  depend¬ 
ing  on  composition.  The  parametric  method  was  used  to  es¬ 
tablish  the  phase  boundaries  between  the  low-temperature 
sphalerite  and  thiogaliate  phases. 

The  parametric  method  takes  advantage  of  the  variation  of 
the  lattice  constants  with  solute  concentration  in  solid  solu¬ 
tions.  It  is  a  classical  technique  which  is  thoroughly  described 
by  Cullity,"  and  requires  no  further  discussion  except  to  note 
that  NaCl  powder  (99.98%)*  was  mixed  with  the  reacted 
samples  for  use  as  an  internal  standard  in  the  X-ray  diffrac¬ 
tion  analysis.  Ni-filtered  Cu Ka  X-rays  were  employed,  with 
the  diffractometer  operating  in  a  digital  step-scanning  mode 
(0.01°  for  10  or  20  s). 

The  basis  of  the  disappearing-phase  method  is  that  the 
fraction  of  a  particular  phase  present  in  the  alloy  is  propor¬ 
tional  to  the  integrated  intensity,  /,  of  one  of  its  X-ray  peaks. 
Klug  and  Alexander"  have  shown  that  in  a  two-phase  mixture 
of,  say,  the  wurtzite  (w)  and  thiogaliate  ( t )  phases,  the  X-ray 
intensity  of  wurtzite  can  be  expressed  as 


_ 

LP  XJq,*  _  at)  +  at  X*fxw  +  A ',ii, 


(1) 


where  /„  is  the  intensity  of  a  specific  wurtzite  peak  in  the 
two-phase  mixture,  Iwp  is  the  intensity  of  the  same  peak,  but 
when  the  material  is  100%  wurtzite  containing  its  equilibrium 
concentration  of  Ga;5-,  (this  is  not  pure  wurtzite),  the  gi’s  are 
mass  absorption  coefficients  of  the  two  phases  at  their  equi¬ 
librium  concentrations  of  GajSi.  and  the  Af’s  are  the  fractions 
by  weight  of  the  two  phases. 

Equation  (1)  is  not  in  an  especially  useful  form  bccaue  /*,,  is 
difficult  to  measure  (measurement  of  lKp  requires  the  prepara¬ 
tion  of  wurtzite  containing  its  equilibrium  concentration  of 
GaiSi  at  a  particular  temperature,  which  is  precisely  what  we 
are  trying  to  determine).  To  transform  it  into  a  useful  equa¬ 
tion,  we  first  write  it  as 


/.  _  LpX^ 

L  +  1 1  lwPXH/xw  f  i,p. X,fi, 


(2) 
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Table  I.  Solid-State  Reaction  Conditions 


no 

Time  (h) 

700 

17(H) 

750 

550 

800 

14(H) 

850 

910 

900 

385 

which  makes  use  of  the  counterpart  of  Eq.  (It  written  for  the 
thiogaliate  phase  instead.  We  then  apply  the  lever  rule.  i.e. 

If;  -  IK,  IK,  -  IK„ 

Y  _ ' _ -  x  =  _ _ _ _ 

'  *  ik  -  w;  '  ik  -  n;. 


where  If;  is  the  weight  fraction  of  Ga.Si  in  the  sample,  and 
IK,  and  W,  are  the  equilibrium  weight  fractions  of  Ga2S5  in 
the  wurtzite  and  thiogaliate  phases,  respectively.  Substitution 
of  these  into  Eq.  (2)  produces  the  result 


IK  -  IK, 


K  +  l'  IK  -  IK,  +  ~-{W„  -  IKJ 

'-eAU 


(3) 


Since  the  ratio  Iv^.l/Iwfp,.  is  constant  within  the  two-phase 
W  +  T  region  at  any  specific  temperature,  and  IK,  <  IK,  < 
If;,  the  ratio  /„/(/„  +  /,)  clearly  varies  between  0  and  1  as  the 
equilibrium  solvus  compositions  are  approached.  If 
were  equal  to  unity,  /„/(/„  +  /,)  would  be  equal  to  X*  and  its 
variation  with  If;,  would  be  linear.  Such  is  not  generally  the 
case,  however,  which  introduces  curvature  into  such  plots 
even  if  the  diffraction  peaks  chosen  are  nearly  equally  intense 
when  the  relative  proportions  of  the  two  phases  are  approxi¬ 
mately  equal  (i.e.,  when  X*  =  X,  ~  0.5). 

In  the  application  of  Eq.  (3)  to  our  experimental  data,  the 
{331}  sphalerite  peak,  the  {0002},  {2023},  {2130}.  and  {2131} 
wurtzite  peaks,  and  the  {112},  {316},  and  {413}  thiogaliate 
peaks  were  used  for  the  intensity  measurements.  The  inte¬ 
grated  intensities  were  calculated  using  the  pow  computer 
program.1"  This  program  subtracts  the  Ka:  component  and 
calculates  the  integrated  intensities  with  a  relative  erroi  of 
less  than  2%. 


III.  Results  and  Discussion 

(l)  Phase  Fields  in  the  ZnS~ZnGa2S4  Pseudobinary  System 
An  analysis  of  the  diffraction  peaks  identified  the  phases 
present  in  the  sample.  The  results  are  plotted  in  Fig.  I,  where 
the  shaded  regions  indicate  the  positions  of  the  phase 
boundaries  consistent  with  the  data.  Preliminary  examination 
revealed  the  existence  of  three  different  solid-solution  phases, 
namely,  sphalerite  (S),  wurtzite  (W),  and  the  tetragonal 
thiogaliate  (T)  and  three  two-phase  fields  (S  +  W,  S  +  T. 
W  +  T).  It  is  obvious  that  the  wurtzite  solid  solution  must 
undergo  the  eutectoid  transformation  W  — »  S  +  T.  but  the 
eutectoid  temperature  and  composition  cannot  be  determined 
accurately  from  the  information  in  Fig.  1.  All  of  these  general 
features  were  noted  by  Malcvskii.4  although  his  diagram  is 
inconsistent  with  the  data  in  Fig.  1. 


Fig.  I.  Phase  map  of  the  ZnS  (ia.S,  system  to 
50  mol';  The  shaded  regions  show  l  he  most  likely  lo¬ 
cations  of  the  phase  boundaries.  The  open  circles  ,n 
dicate  the  presence  of  sphalerite  (S).  the  open  squares 
the  presence  of  wurt/ite  (Wl,  and  the  filled  circles 
the  presence  ol  thiogaliate  (T). 
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(2)  Quantitative  Evaluation  of  the  Phase  Boundaries 

Figure  2  shows  examples  of  the  use  of  the  disappearing- 
phase  method  used  to  determine  the  S/W  and  W/T  solvi  at 
two  different  temperatures.  The  intensity  ratios  /,/(/,  +  /., ' 
and  /»/(/»  +  /,)  are  plotted  vs  wt%  Ga:Si  for  the  S  f  W  and 
W  +  T  phase  fields,  respectively.  The  use  of  weight  percent 
rather  than  mole  percent  was  dictated  by  the  fact  that  Eq.  (3) 
is  derived  expressly  for  these  concentration  units.  The  data 
were  fitted  by  a  fourth-order  polynomial  for  purposes  of  ex¬ 
trapolation  to  0%  and  100%  of  the  phases  present,  and  the 
intersections  were  taken  as  the  equilibrium  compositions. 
The  use  of  a  high-order  polynomial  obviated  the  need  to  im¬ 
pose  any  assumptions  on  the  behavior  of  the  intensity  ratios 
with  composition,  and  illustrates  the  nonlinear  nature  of 
Eq.  (3). 

At  800°C  and  below,  the  two-phase  region  consists  only  of 
sphalerite  and  thiogallate,  and,  as  already  noted,  the  paramet¬ 
ric  method  is  more  accurate  in  this  temperature  range.  Ex¬ 
amples  of  its  use  are  shown  in  Fig.  3.  For  the  ZnS-rich  end, 
the  lattice  constant  of  sphalerite,  a,  is  plotted  vs  mol%  Ga2S, 
using  the  angular  position  of  the  {331}  sphalerite  peak  to  cal¬ 
culate  a.  The  solvus  compositions  at  700°  and  750°C  are  lo¬ 
cated  at  5.1  and  6  mol%  GajS.v  At  the  ZnGa2S4-rich  end,  the 
positions  of  the  {400},  {316},  and  {413}  peaks  were  used  to 
calculate  a,  as  a  function  of  the  mol%  Ga2Si.  The  boundary 
between  the  one-  and  two-phase  regions  is  44%  Ga^S3  at 
700  C  and  41%  at  800°C. 

The  ZnS-Ga2S3  phase  diagram  thus  obtained  from  a  com¬ 
bination  of  the  two  methods  is  shown  in  Fig.  4.  Extrapolation 
of  the  phase  boundaries  between  the  one-  and  two-phase 
fields  was  used  to  determine  the  eutectoid  composition  and 
temperature,  which  were  assumed  to  be  the  point  of  inter¬ 
section  of  the  high-temperature  data.  Since  the  data  at  high 
temperatures  were  limited,  the  extrapolation  produces  results 
which  are  accurate  only  to  within  about  1  mol%.  Neverthe- 


(a) 


Fig.  2.  Illustrating  the  use  of  the  disappearing-phase 
method  to  determine  the  boundaries  between 
(a)  wurtzite  and  sphalerite  at  XS0  C  and  (b)  wurt/.ite 
and  thiogallate  at  900C'. 


Fig.  3.  Illustrating  the  use  of  the  parametric  method  to 
determine  (a)  the  S/(S  +  T)  solvus  at  700°  and  75(TC  and 
(b)  the  T/(S  +  T)  solvus  at  700°  and  800"C.  In  (b).  a 
refers  to  a,. 

less,  the  coordinates  of  the  intersection  point,  namely,  16% 
Ga2Si  and  818°C,  represent  the  best  estimate  of  the  eutec¬ 
toid  composition  and  temperature  currently  available.  A  simi¬ 
lar  extrapolation  applied  to  the  phase  boundaries  at  the 
sphalerite-rich  and  thiogallate-rich  ends  of  the  phase  diagram 
places  the  compositions  of  the  reactant  phases  at  —9  and 
—40  mol%  Ga:Sj,  respectively.  The  phase  boundaries  indi¬ 
cated  by  the  shaded  curves  in  Fig.  1  arc  consistent  with  those 
determined  quantitatively.  It  is  difficult  to  determine  how  the 
errors  in  all  the  measurements  propagate  to  the  compositions 
and  temperatures  cited.  We  estimate  that  the  eutectoid  tem¬ 
perature  is  accurate  to  ±5°C  and  the  compositions  accurate 
to  ±1  mol%  Ga2Si. 

On  comparing  our  phase  diagram  with  that  of  Malevskii,4 
we  note  that  the  accuracy  of  his  diagram  is  limited  by  the 
much  shorter  reaction  times  used  (400  h  at  650°C;  cf.  more 


Fig.  4.  ZnS  (ia-Si  phase  diagram  to  sc  ,  oh  i .  from 
700  to  000  ( 
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than  161X1  h  a!  7(X)~C  in  our  ease),  and  the  much  larger  inter¬ 
vals  in  composition  employed,  especially  near  the  eutectoid 
composition  (roughly  7  molrf).  These  two  factors  are  un¬ 
doubtedly  responsible  for  the  discrepancies  between  the  two 
phase  diagrams. 

We  note  that  the  suiuHlify  of  ZnS  increases  with  increasing 
temperature.  Consequently,  both  eutectoid  decomposition 
and  precipitation  hardening  can  be  potentially  utilized  to  ma¬ 
nipulate  the  mechanical  properties  of  ZnS-Ga:S1  ceramics, 
as  discussed  by  Zhang  et  al.2 


IV.  Summary 

The  ZnS  -Ga  -S,  phase  diagram,  from  7(X)°  to  9(XTC  and 
to  50  moK'f  Ga:S-,,  was  determined  using  X-ray  diffraction 
techniques  from  samples  prepared  via  the  solid-state  reaction 
of  powders.  The  wurtzite  (W),  sphalerite  (S),  and  tetragonal 
thiogallate  (T)  phases  all  exist  over  a  considerable  range 
of  stoichiometry.  There  is  a  eutectoid  reaction  at  818  ±  5°C 
in  which  W  (16  ±  1  mol%)  decomposes  to  S  (9  ±  1%)  + 
T  (40  ±  lr7).  The  solubility  of  Ga^Sj  in  ZnS  increases  with 
increasing  temperature  below  the  eutectoid  temperature. 
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ABSTRACT: 

The  ZnS-rich  region  of  the  ZnS-CdS  phase  diagram,  up  to  18  mol  %  CdS,  was 
studied  over  the  temperature  range  700  to  950  °C,  and  the  boundaries  describing  the 
equilibria  between  the  sphalerite  and  wurtzite  solid  solutions  were  determined  X-ray 
diffraction  analysis  was  the  principal  technique  used,  the  phase  boundaries  being  de¬ 
termined  by  the  lattice-parametric  and  disappearing  phase  methods.  The  re  ults  indi¬ 
cate  that  the  solubility  of  CdS  in  ZnS  is  extensive  in  both  the  sphalerite  and  wurtzite 
phases,  with  CdS  stabilizing  the  wurtzite  phase. 

MATERIALS  INDEX:  zinc,  sulfides,  cadmium 


Introduction 


The  II-VI  compound  ZnS  is  technologically  important  due  to  its  infrared  transmitting  and 
luminescent  properties.  Although  many  of  its  optical  properties  are  adequate  for  various  applica¬ 
tions,  there  is  considerable  interest  in  improving  its  mechanical  behavior.  One  method  for  harden¬ 
ing  and  toughening  ceramic  materials  is  by  the  incorporation  or  precipitation  of  a  second  phase  in 
the  matrix.  This  can  be  accomplished  by  a  variety  of  methods,  including  thermal  or  thermo-me¬ 
chanical  processing  to  control  solid-state  reactions,  which  is  the  route  we  are  pursuing.  To  utilize 
this  approach  a  know. edge  of  the  phase  diagrams  of  the  systems  of  potential  interest  is  necessary. 

ZnS  is  known  to  exist  in  two  polymorphic  modifications,  sphalerite  (cubic,  zinc-blende) 
and  wurtzite  (hexagonal).  The  cubic  form  is  stable  at  low  temperatures  and  transforms  to  wurtzite 
at  1020  °C  (1).  The  sphalerite-wurtzite  (S-W)  transition  is  affected  by  the  addition  of  solutes,  but 
knowledge  of  the  phase  equilibria  is  sparse.  We  selected  the  ZnS-CdS  system  for  study  because 
the  optical  properties  of  CdS  are  such  that  ZnS-CdS  solid  solutions  or  two-phase  mixtures  are 
likely  to  retain  the  excellent  infrared  transmission  characteristics  of  pure  ZnS. 

CdS  forms  a  complete  series  of  solid  solutions  with  wurtzite  at  high  temperatures  (T  >  975 
°C)  (2),  but  the  only  work  on  the  phase  equilibria  at  lower  temperatures  is  a  study  conducted  under 
hydrothermal  conditions  over  the  temperature  range  250-450  °C  (3).  Those  conditions  and  tem¬ 
peratures,  however,  are  not  practical  for  the  thermal  processing  treatments  we  believe  are  required 
to  produce  strengthening  and  toughening.  In  this  investigation,  the  phase  boundaries  between  the 
ZnS-rich  sphalerite  and  wurtzite  solid  solutions  were  determined  over  the  temperature  range  700- 
950  °C. 


1667 


1068 


Y< .  W.  CHEN,  et  al. 


Vol.  23,  No.  11 


Experimental 


Samples  were  prepared  from  weighed  amounts  of  cubic  ZnS  (99.9  9c,  Aesar)  and  hexago¬ 
nal  CdS  (99.9  +  %.  Aesar)  powder  and  mixed.  The  typical  sample  weighed  at  least  1  _,  and  alios., 
containing  up  to  18  mol  9c  CdS  were  investigated.  Approximately  0.2  p  of  each  mixture  w'as 
transfened  to  a  silka  ampoule,  evacuated  to  2  x  10'2  torr,  and  heated  to  120  °C  for  10  ‘o  15  min¬ 
utes  to  remove  adsorbed  water.  The  ampoule  was  then  sealed  and  isothermallv  annealed  at  the  de¬ 
sired  temperature  for  times  ranging  from  1  to  4  weeks.  X-ray  diffraction  analysis  using  mono¬ 
chromatic  Ni-filtered  CuKa  radiation  was  employed  to  identify  the  phases  present  and  to  determine 
their  lattice  constants.  The  phase  boundaries  were  evaluated  by  using  either  the  disappearing  phase 
method  or  the  parametric  method  M).  The  X-ray  intensities  were  recorded  digitally  in  a  step-scan¬ 
ning  mode  at  intervals  of  0.01°  for  1  s  or  0.02°  for  10  s.  Eithe1'  NaCl  or  KC1  powders  were  used 
as  a  standard,  depending  upon  the  amount  of  peak  overlap  wun  the  samples,  to  obviate  the  need 
for  precise  alignment  of  the  diffractometer. 

R'  ults  and  Discussion 

The  phases  detected  by  analysis  of  the  peaks  in  the  X-ray  diffraction  patterns  in  a  large 
number  of  samples  heat-treated  at  the  four  temperatures  used  in  this  study  are  presented  in  Fig.  1 . 
This  is  a  map  describing  the  phase  equilibria,  in  which  the  compositions  of  the  phase  boundaries 
are  only  approximately  determined.  The  compos:aons  w-ere  then  determined  more  precisely  using 
the  X-ray  methods  mentioned  above. 
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FIG.  1 

Identification  of  the  phases  observed  on  heating  mixtures  of 
ZnS  and  CdS  powders  for  various  times  at  the  temperatures  indicated. 


An  example  of  the  determination  of  a  phase  boundary  using  the  parametric  mctl  ,»d  is 
shown  in  Fig.  2.  In  a  single-phase  region  the  lattice  constant  of  the  solid  soluuon  will  vary  with 
the  concentration  of  solute  if  the  atomic  or  ionic  radii  of  the  solvent  and  solute  differ.  Generally, 
this  change  is  approximately  linear  .rnd  should  be  independent  of  the  thermal  history  of  the  sample 
so  long  as  equilibrium  is  reached.  The  solid  line  in  Fir.  2  represents  the  variation  of  the  lattice 
constant  with  composition  of  the  sphalerite  solid  solutioi  in  the  single-).  ^a,e  region.  At  ihe  higher 


Vol .  23,  No.  II 


ZnS-CdS  SYSTEM 


1669 


reactioii  temperatures  (900  and  950  °C)  equilibrium  was  established  relatively  rapidly.  The  lattice 
constant  of  spha'crn.  is,  of  course,  invariant  in  the  two-p  iase  S  +  W  field,  but  the  data  are  limited 
by  the  fact  that  the  phase  field  is  narrow.  Hence,  the  number  of  data  points  at  these  two  tempera¬ 
tures  is  small.  As  is  evident  in  Fig.  2.  the  lattice  constants  for  the  two-phase  mixtures  in  samples 
reacted  at  the  lower  temperatures  (700  and  800  °C)  continued  to  increase  with  increasing  concen¬ 
tration  of  CdS,  but  not  as  rapidly  as  in  the  single-phase  solid  solution.  This  is  probably  due  to  the 
incomplete  establishment  of  equilibrium.  Consequently,  the  accuracy  of  the  sphalerite  solvus  dc 
termined  from  the  parametric  method  is  not  as  valid  as  those  at  tKK)  and  950°C  due  to  the  variance 
of  the  lattice  constant  values  in  the  two-phase  region.  For  this  reason  we  relied  on  the  disappearing 
phase  methc-l  to  determine  the  phase  boundaries  at  the  lower  temperatures. 
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FIG.  2 

The  variation  of  the  lattice  constant  of  sphalerite,  a().  as  a  function  of  the 
C  IS  content  of  the  mixtures  prepared.  Tht  points  indicated  by  x  apply  to  the 
single  phase  solid  solution  of  CdS  in  ZnS. 


The  parametric  method  was  also  used  to  determine  the  (S  +  W)/W  phase  boundary  com¬ 
pose  on  s  at  900  and  950  °C.  The  c  and  a  lattice  constants  both  vary  lirearlv  with  the  CdS  concen¬ 
tration  in  the  single-phase  wurtzite  region,  and  the  intersection  between  the  horizontal  and  sloped 
lines  is  the  solvus  as  show  n  in  Fig.  v  The  compositions  of  the  wurtzite  phase  boundaries  deter¬ 
mined  from  both  c  and  a  parameters  are  consistent  (-4  mol  %  CdS  at  950  °C  and  5.5  to  6  mol  G 
CdS  at  900  °C).  The  small  number  of  data  points  determining  the  horizontal  lines  is  due  to  the 
narrow  compositional  range  of  the  two-phase  field  at  these  temperatures. 

The  disappearing  phase  method  proved  to  be  more  reliable  for  determining  the  Dhase 
boundaries  at  the  lower  temperatures.  The  method  is  based  upon  the  variation  of  the  relative 
amounts  of  S  and  W  along  a  tie-line  in  the  two-phase  field  as  the  overall  composition  is  changed 
(4).  For  example,  the  weight  fraction  of  W  should  decrease  to  0  as  the  compt  rition  appir  aches  the 
S  solvus.  and  since  the  intensity  of  a  particular  wurtzite  peak,  Iw.  is  propon  >nal  to  the  weight  G 
W  in  the  two-phase  mixture,  the  ratio  Iw/Is  should  also  approach  a  value  of  zero  a.,  the  sphalerite 
solvus  is  reached  (Is  is  the  intensity  of  a  sphalerite  peak).  The  solvus  composition  is  determined 
from  an  extrapolation  of  the  variation  of  Iw/Is  vs.  weight  %  CdS  to  the  zero  ordinate  value. 
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FIG.  3 

Results  obtained  applying  the  parametric  method;  (a)  illustrates  the  variation  of  the  lattice  constant 
a,  while  (b)  illustrates  the  variation  of  the  lattice  constant  c  with  the  concentration  of  CdS. 


The  integrated  intensities  of  the  (31 1 )  sphalerite  and  the  (2131 )  wurtzite  peaks  were  cho¬ 
sen  for  this  analysis.  Is  and  Iw  were  computed  using  a  numerical  analysis  program,  subtracting  the 
contribution  of  CuKa2  radiation.  Depending  on  the  system  being  investigated,  the  relative  ratio  of 
the  intensities  can  vary  linearly  or  non-linearly  with  composition  (4).  From  a  cursory  examination 
of  the  data,  we  decided  to  try  a  linear  relationship  and  compare  the  results  with  the  data  in  Fig.l. 
The  data  are  shown  in  Fig.  4,  where  it  is  evident  that  the  assumed  linear  behavior  is  better  obeyed 
by  the  data  in  Fig.  4a  than  in  Fig.  4b. 

From  a  combination  of  the  two  analytical  methods  used,  the  partial  phase  diagram  for  the 
ZnS-CdS  system  over  the  temperature  range  of  this  study  is  shown  in  Fig.  5a.  The  solid  curves 
drawn  in  Fig.  5a  are  shown  in  Fig.  5b,  where  they  are  superimposed  on  the  data  from  Fig.  1;  the 
agreement  is  very  good.  We  are  therefore  confident  that  the  results  obtained  using  the  disappearing 
phase  method  are  accurate,  despite  the  scatter  in  Fig.  4b.  The  disappearing  phase  method  suc¬ 
ceeded  where  the  parametric  method  failed  because  the  integrated  intensities  of  the  broad  peaks  in 
the  two-phase  mixtures  are  more  accurately  measured  than  the  positions  of  the  peak  maxima 

The  partial  ZnS-CdS  phase  diagram  determined  over  the  temperature  range  250  to  450  °C 
under  hydrothermal  conditions  by  Kaneko  et  al.  (3)  exhibits  features  similar  to  those  in  Fig.  5,  es¬ 
pecially  insofar  as  the  solubility  of  CdS  in  ZnS  decreases  with  increasing  temperature.  Kaneko  et 
al.  made  no  attempt  to  determine  the  phase  boundaries,  but  it  is  evideni  that  their  two-phase  equi¬ 
libria  are  shifted  to  higher  mol  %  CdS,  and  the  compositional  range  of  their  two-phase  field  is 
much  greater.  The  two  phase  diagrams  can  be  reconciled  only  if  the  sphalerite  phase  is  stabilized 
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under  hydrothermal  conditions,  but  there  is  no  convincing  evidence  to  suggest  that  this  stabilization 
process  occurs.  Cherin  et  al.  (5)  found  a  two-phase  mixture  of  sphalerite  and  wurtzite  in  their  ZnS 
single  crystals  containing  5  mol  %  CdS.  The  crystals  were  grown  in  a  temperature  gradient  of 
925  to  950  °C  using  the  iodine  vapor  transport  method.  Their  result  is  entirely  consistent  with  our 
phase  diagram. 


5  6  7  8  9  10  11121314 


wt  %  CdS 


FIG.  4 

Results  obtained  using  the  disappearing  phase  method:  (a)  The  ratio  of  the  intensities  of  the 
wurtzite  (2151  |  peak  and  the  sphalcrm  1 3 !  1 }  peak,  lw/ls,  as  a  tunction  of  the  concentration  of 
CdS  in  weight  91.  In  (hi  the  ratio  of  the  intensities  is  inverted. 


The  use  of  differential  thermal  analysis  for  determination  of  the  phase  boundaries  was  at¬ 
tempted  in  several  trial  runs  using  pure  ZnS.  Since  the  enthalpy  of  the  sphalerite  to  wurtzite  trans¬ 
formation  is  3.19  kcal/mol  (6)  in  pure  ZnS.  the  transition  was  not  detected  on  heating  and  cooling 
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through  1020  °C.  Thus,  the  small  heat  evolution  during  the  transformation  rendered  thermal 
analysis  ineffective  for  determining  the  phase  diagram  in  this  system. 


1100 
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mol  %  CdS 

FIG.  5 

(a)  Partial  phase  diagram  for  the  ZnS-CdS  system  over  the  temperature  range  700  to  1 100  °C.  A 
comparison  with  the  phases  observed  over  the  range  of  conditions  in  Fig.  1  is  shown  in  (b). 

■Summary 

The  partial  phase  diagram  for  the  ZnS-CdS  system  was  determined  over  the  temperature 
range  700  to  950  °C.  Cubic  and  hexagonal  solid  solutions  of  ZnS  were  formed,  and  the  bound¬ 
aries  of  the  two  phase  region  were  determined  by  the  use  of  the  parametric  and  disappearing  phase 
methods. 
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ABSTRACT 

A  series  of  ZnS-rich  wurtzite  solid  solutions  containing  8,  12  and  16  mol%  Ga2S3has  been  prepared.  The 
influences  of  the  fabrication  conditions  and  microstructure  on  the  mechanical  properties  were  examined.  It  has  been 
demonstrated  that  Zn|_xGa2XS[  +2xis  tougher  and  harder  than  pure  ZnS  due  to  a  solid-solution  strengthening  effect. 
The  hardness  of  the  solid  solutions  increases  with  increasing  concentration  of  Ga2S3and  decreasing  grain  size,  the 
grain-size  dependence  being  of  the  Hall-Petch  type.  Finer-grained  solid  solutions  with  moderate  concentrations  of 
Ga2S3have  the  highest  fracture  toughness.  Whereas  the  hardness  increases  monotonically  with  increasing  mole 
fraction  of  Ga2S3,the  fracture  toughness  is  a  maximum  at  around  12  %  Ga2S3. 

1.  INTRODUCTION 

ZnS  is  used  as  an  infrared  window  material  for  applications  in  the  8-12  gm  range  because  of  its  adequate  optical 
transmission  and  mechanical  behavior.  However,  improvements  in  its  fracture  toughness,  thermal  shock  resistance 
as  well  as  erosion  resistance  are  highly  desirable  for  missile  dome  applications.  To  obtain  useful  optical  properties 
and  improved  mechanical  behavior  requires  a  microstructure  which  is  inhomogeneous  with  regard  to  deformation 
and  fracture,  but  homogeneous  for  the  wavelength  of  interest.  Since  multiphase  materials  are  inherently  undesirable 
for  optical  transmission,  attempts  were  made  to  produce  single-phase  solid  solutions  of  ZnS  (or  the  solid  solution 
plus  very  fine  precipitates)  which  are  not  only  infrared  transmitting  in  the  8  to  12  gm  range  but  also  possess 
improved  mechanical  properties.  Solid-solution  strengthening  has  been  found  in  some  ceramic  systems1,2.  In  this 
study,  ZnS-Ga2S3  solid  solutions  were  synthesized  to  increase  the  mechanical  strength  of  ZnS.  Ga2S3 was  chosen 
as  a  second  phase  component  for  reasons  discussed  by  Zhang  et  al2. 

2.  EXPERIMENTAL 

ZnS-x  mol  %  Ga2S3  solid  solutions  (where  x  is  8,  12  and  16  respectively)  were  synthesized  by  sintering  a 
physical  mixture  of  powders  of  ZnS  (99.99%,  Aldrich)  and  Ga2S3  (99.99%,  Alfa).  The  solid-state  reaction  was 
completed  after  18  hours  at  925 °C  in  H2S  atmosphere.  The  powders  were  identified  as  single-phase  wurtzite 
(hexagonal,  B4  structure)  by  X-ray  diffraction  analyses.  A  vibromill  (8000  Mixer/Mill,  SPEX)  was  used  to  reduce 
the  site  of  the  starting  solid  solution  powders  to  about  3  /im.  The  powder  was  cold-pressed  in  vacuum  at  207  MPa 
and  was  subsequently  hot-pressed  in  a  molybdenum  die.  Table  1  lists  the  hot-pressing  conditions  (Tp,  t  and  p  are 
the  hot-pressing  temperature,  holding  time  and  pressure  respectively)  for  all  alloys  and  pure  ZnS  and  Fig.  1  shows 
a  schematic  hot  pressing  schedule.  The  sample  was  heated  at  300°C  for  10  minutes  to  eliminate  water  adsorbed 
on  the  powders,  then  continuously  heated  at  a  constant  rate  of  20°C/min  to  900°C,  after  which  the  pressure  (138 
MPa)  was  applied.  Temperature  and  pressure  were  held  for  about  20  to  40  minutes  and  then  usually  the  pressure 
was  released  and  sample  slowly  cooled  to  room  temperature.  The  densities  of  the  hot-pressed  pellets  were  measured 
using  a  liquid  immersion  technique  based  on  the  ASTM  standard4. 

A  series  of  hot-pressed  pellets  was  sealed  into  quartz  ampoules  and  solid-solution  treated  in  the  wurtzite  phase 
field  for  various  lengths  of  time.  The  heat-treatment  temperatures  were  selected  according  to  the  partial  ZnS-Ga2S3 
phase  diagram  we  established  earlier2.  The  purpose  of  the  treatments  is  twofold:  (a)  to  ensure  that  the  specimen 
is  completely  transformed  to  wurtzite;  (b)  to  grow  grains  of  different  sizes.  Pure  ZnS  was  annealed  at  1080°C  and 
all  three  solid  solutions  at  940°C.  The  ampoules  Were  water-quenched  to  room  temperature  immediately  after  all 
reaction  treatments.  Microprobe  analysis  and  EDS  were  performed  on  ZnS  and  all  the  solid  solutions,  and  the 
results  confirmed  the  expected  compositions. 
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The  hardness  and  fracture  toughness  were  evaluated  using  the  Vickers  indentation  technique,  the  hardness, 
H,  in  GPa,  is  obtained  using  the  formula 

H  =  18.169—  (1) 

a1 

where  P  is  the  applied  load  in  grams  and  a  is  the  diagonal  length  of  the  Vickers  indentation  in  fi m. 

Table  1.  Hot  Pressing  Parameters. 
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Kig.  I.  Hot  pressing  program 


Fracture  toughness  measurement  by  the  Vickers  indentation  technique  is  based  on  indentation  fracture 
mechanics'’.  Because  of  its  uniquely  simple  and  economic  nature,  this  technique  has  been  widely  used  for  measuring 
fracture  toughness,  especially  for  small  specimens.  The  fracture-toughness  data  obtained  by  the  indentation 
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technique  are  accurate  to  within  either  10%  if  Young’s  modulus  is  known  or  30%  if  Young’s  modulus  is  not 
known6.  The  fracture  toughness  for  a  median/radial  crack  system,  /qc,  can  be  expressed  as7 


■  «|)'B 


P 

cm 


(2) 


where  6  is  a  material-independent  constant  (-0.016),  E  is  Young’s  modulus  and  C  is  the  crack  length. 

Specimens  approximately  one  quarter  of  the  12.4  x  2  mm  discs  were  first  cold-mounted  in  epoxy  and  then 
ground  and  polished.  Final  polishing  was  done  on  a  standard  metallographic  wheel  using  0.05  fim  alumina 
polishing  powder.  Considerable  care  was  taken  to  ensure  smooth  and  residual-stress  free  surfaces  to  measure 
hardness  and  toughness.  A  Microhardness  Tester  (MicroMet  11,  Buehler)  was  used  to  measure  H  for  at  least  three 
different  loads  (typically  200,  500  and  1000  g)  and  more  than  18  indentations  for  each  specimen  were  taken  using 
a  loading  time  of  15  s.  The  indentation  impression  size  and  Vickers  indentation  crack  length  were  measured 
immediately  after  unloading.  It  was  observed  that  cracks  propagate  in  an  intragranular  mode  for  all  the  specimens. 
Young’s  modulus,  E,  for  all  the  specimens  was  measured  using  a  Nanoindenter,  yielding  an  average  value  of  90 
GPa.  This  value  was  used  in  the  calculations  of  ^according  to  equation  (2). 

3.  RESULTS  AND  DISCUSSION 

Since  the  highest  pressure  is  limited  by  the  strength  of  the  die  material,  it  was  necessary  to  explore  different 
conditions  of  Tp,  p  and  t  for  hot  pressing.  The  optimum  hot-pressing  condition  for  the  ZnS-GaiS^  solid  solutions 
is  Tp  =  900°C,  p  =  124  to  138  MPa  and  t  =  20  to  40  minutes.  The  conditions  used  to  produce  specimens 
exceeding  99.0%  of  their  theoretical  densities  are  listed  in  Table  1. 
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Fig.  2.  Hardness  vs.  the  reciprocal  square 
root  of  the  grain  size  for  ZnS  and 
the  solid  solutions 


%  00 


Mol  »  Ga-jS-, 

Fig.  3.  Hardness  vs.  concentration  of  Ga-S;, 
The  curve  A,  B  and  C  were  obtained 
from  the  curves  fitted  to  the  data 
in  Fig.  2. 


The  hardness  values  of  the  specimens  as  a  function  of  grain  size,  d,  are  listed  in  Table  2,  and  plotted  as  a 
function  of  d  in  Fig.  2.  It  can  been  seen  that  hardness  increases  as  d  decreases  for  all  three  solid  solutions  and 
pure  ZnS.  This  indicates  that  smaller  grains  produce  higher  resistance  to  dislocation  motion,  possibly  due  to 
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dislocation  pile-ups  at  grain  boundaries.  Since  the  hardness  of  a  ceramic  material  is  aboui  three  times  its  yield 
stress*,  the  dependence  of  H  on  d  can  be  expressed  using  the  Hall-Petch  equation: 

H  =  H0  +  kd'in  W 


where  Hq  and  k  are  constants.  The  present  experimental  data  are  in  good  agreement  with  equation  (3). 

Table  2.  Hardness  and  toughness  of  pure  ZnS  and  the  ZnS-Ga2S3 
alloys  as  a  function  of  grain  size. 


1 

Materials 

Grain 

Diameter 

(nm) 

Hardness 

(GPa) 

Toughness 

(MN/m3'2) 

ZnS 

5 

2.62  +  0.09 

0.50  ±  0.10 

100 

1.44  +  0.23 

0.40  ±  0.09 

ZnS-8%  Ga->S3 

6 

3.59  +  0.14 

0.72  +  0.09 

16 

3.11  +0.16 

0.68  +  0.10 

18 

3.09  ±  0.20 

0.66  ±  0.15 

ZnS- 12%  Ga2S3 

6 

4.27  +  0.27 

0.77  +  0.1S 

14 

4.01  ±  0.27 

0.70  ±  0.12 

ZnS- 16%  G22Sj 

7 

4.25  +  0.17 

0.62  +  0.09 

16 

4.09  +  0.22 

0.54  +  0.10 

20 

3.95  +  0.17 

0.53  +  0.09 

The  hardness  as  a  function  of  the  amount  of  Ga2S3  in  solid  solution  is  plotted  in  Fig.  3.  The  data  were 
extracted  from  Fig.  2,  with  curve  D  obtained  by  extrapolating  to  H  =  H0  at  du2  =  0.  It  is  clear  from  Fig.  3  that 
the  hardness  of  the  solid  solutions  increases  monotonically  with  Ga2S3 content.  The  alloy  with  16  %  Ga^was 
about  twice  as  hard  as  pure  ZnS  in  the  grain  size  range  studied.  The  increase  in  the  hardness  for  the  solid  solutions 
results  from  the  interactions  between  dislocations  and  defects  induced  by  the  Ga2S3 dissolved  in  the  crystals. 

The  resistance  to  dislocation  motion  that  constitutes  solid-solution  strengthening  can  come  from  one  or  more 
of  the  following  factors.  First,  a  stress  field  will  be  induced  by  the  relative  size  difference  between  Ga  and  Zn  ions; 
the  actual  difference  is  about  25  % ,  the  sizes  of  Zn  and  Ga  ions  being  0.083  and  0.062  nm  respectively.  In  addition, 
the  substitution  of  two  Ga+3ions  for  three  Zn  +  2will  create  a  vacancy  to  preserve  electrical  neutrality,  and  the 
hardening  can  be  increased  by  an  associated  vacancy -Ga+ 3  defect.  Additionally,  a  difference  in  the  effective  shear 
modulus  between  Ga  and  Zn  could  be  an  important  factor,  and  an  electrical  interaction  between  dislocations  and 
Ga+3ions  or  vacancy-Ga+3complex  could  impede  dislocation  movement.  Finally,  the  spatial  distribution  of  the 
solute  ions  might  not  be  completely  random,  and  either  short-range  ordering  or  clustering  would  produce  a 
hardening  effect.  Many  of  these  mechanisms  are  known  to  be  effective  in  inhibiting  the  motion  of  dislocations  in 
metallic  alloys  and  can  be  expected  to  operate  in  ceramic  alloys  as  well.  Identifying  the  specific  mechanism(s) 
requires  additional  research,  however. 

Fracture  toughness  for  all  the  specimens  is  plotted  versus  rZ12in  Fig.  4.  It  is  readily  seen  that  all  the  solid 
solutions  are  tougher  than  ZnS,  with  ZnS- 12  %  Ga2S2  having  the  highest  toughness  values.  Figure  4  also  shows 
that  Kjc  is  approximately  proportional  to  d 1/2  for  both  pure  ZnS  and  the  solid  solutions  in  the  grain  size  range 
covered.  The  cracks  emanating  from  the  microhardness  indentations  are  shown  in  Fig.  5.  where  it  is  evident  that 
the  crack  patterns  in  the  solid  solutions  are  more  irregular  than  those  of  pure  ZnS,  the  cracks  often  curving, 
branching  and  deviating  from  the  plane  intersecting  the  comers  of  the  indentation.  The  micrographs  in  Fig.  5 
indicate  that  propagating  cracks  are  more  easily  deflected  in  the  solid  solutions  than  in  pure  ZnS. 
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Of  all  toughening  mechanisms  discussed  by 
Swain9,  transformation  toughening  is  unlikely  to  be 
of  any  importance  to  this  system  because  there  is  only 
a  0.2  %  volume  change  accompanying  the  transforma¬ 
tion  from  wurtzite  (W)  to  sphalerite  (S)  (the  cubic 
form  of  the  ZnS  solid  solutions,  which  has  a  smaller 
molecular  volume),  though  it  is  possible  to  induce  the 
W-S  transformation  by  an  applied  stress10.  On  the 
other  hand,  microcrack  toughening  is  one  of  the  most 
likely  mechanisms  to  cause  crack  deflection  in  the 
ZnS  solid  solutions.  Microcrack  toughening  has  been 
found  to  be  an  effective  toughening  mechanism  in 
many  ceramic  materials.  In  non-cubic  ceramics 
produced  by  hot  pressing,  microcracks  develop  on 
cooling  from  the  processing  temperature.  According 
to  Rice  et  al. 1 1  microcrack  toughening  by  this  mecha¬ 
nism  is  grain-size  dependent,  the  number  of  microcra¬ 
cks  per  unit  volume  increasing  with  increasing  grain 
size  while  the  elastic  energy  absorbed  per  microcrack 
decreases  with  increasing  grain  size.  This  model 
therefore  predicts  an  optimum  value  of  the  grain  size 
which  maximizes  K[c. 

Lewis  et  al.12  have  shown  that  K Ic  for  CVD- 
grown  cubic  ZnS  behaves  similarly  in  that  Ki{.  is 
maximized  at  some  value  of  d.  They  suggest  that  the 
grain-size  dependence  of  K\c  is  due  to  crack  nucle- 
ation  through  dislocation  pile-ups  at  grain  boundaries, 
and  the  decrease  in  K\c  for  very  small  grain  sizes  is 
attributed  to  void  formation  in  a  spherical  zone 
beneath  the  Vickers  indentation  contact  area.  Al¬ 
though  this  is  not  a  microcrack  toughening  mecha¬ 
nism,  the  dependence  of  K (c  on  d  can  nevertheless  be 
similar  for  cubic  and  non-cubic  materials. 

Whereas  a  maximum  value  of  /qcas  a  function  of 
d  is  not  observed  in  Fig.  4,  the  mechanism  of  Rice  et 
al.  could  nevertheless  be  viable  in  hexagonal  ZnS. 
When  the  non-cubic  polycrystalline  ZnS  solid  solu¬ 
tions  are  cooled  from  the  hot-pressing  temperature  to 
room  temperature,  intercrystalline  microstresses  can 
be  expected  to  develop  as  a  result  of  the  thermal - 
expansion  anisotropy  of  the  grains.  The  mechanism 
proposed  by  Lewis  et  al.12  could  also  contribute  to 
the  increasing  toughness,  however,  because  there  is 
no  reason  why  plasticity  (as  opposed  to  microcrack) 
toughening  cannot  play  a  role  in  augmenting  the 
toughness  of  non-cubic  materials.  At  this  stage  of 
the  research  there  is  no  way  of  evaluating  the  relative 
importance  of  the  microcrack  vs.  plasticity  toughen¬ 
ing  mechanisms. 
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Fig.  4.  Fracture  toughness  vs.  the  reciprocal 
square  root  of  the  grain  size  for  ZnS 
and  the  solid  solutions 


Fig.  5.  Vickers  hardness  indentations  under  a  load  of 
200  g  for  (a)  pure  ZnS  and  (b)  for  ZnS- 16% 
Ga2S3 
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It  can  he  observed  in  Fig.  4  that  the  effect  of  the  Ga2Sj  concentration  on  the  toughness  of  the  solid  solutions 
is  not  a  monotonic  function.  Of  the  three  ceramic  alloys,  ZnS-12%  G 8283  has  the  highest  toughness  values, 
implying  that  there  is  an  optimum  Ga2S3Concentration  for  fracture  toughness.  The  underlying  physical  reason  for 
an  optimal  concentration  of  GajSj  is  currently  not  understood  and  remains  another  interesting  topic  for  future 
research. 


4.  SUMMARY 

1.  The  incorporation  of  Ga2S3into  ZnS  improves  both  the  hardness  and  toughness  of  ZnS. 

2.  The  hardness  of  the  ZnS-Ga2S3  solid  solutions  is  proportional  to  the  amount  of  Ga2S3in  the  solutions. 

3.  The  hardness  increases  as  the  grain  size  decreases  for  all  three  solid  solutions  and  pure  ZnS.  It  is  consistent 
with  the  Hall-Petch  relation:  H  =  Hq+  kd~v~. 

4.  All  the  solid  solutions  are  tougher  than  ZnS.  A  finer  grain  size  results  in  higher  toughness.  However,  the 
dependence  of  fracture  toughness  of  the  solutions  on  the  solute  content  is  not  a  monotonically  increasing 
function.  It  seems  that  there  is  an  optimum  composition  for  the  fracture  toughness. 
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ABSTRACT 


Novel  apparatus  and  methodology  have  been  developed  for  measur¬ 
ing  the  fracture  toughness  of  ceramics  using  small  disk-shaped  samples  3 
mm  in  diameter.  The  method  involves  the  Vickers  indentation  of  speci¬ 
mens  ranging  in  thickness  from  300  to  700  pm,  and  testing  them  in  a  ring- 
on-ring  bending  mode.  Experiments  on  CVD-grown  ZnS  have  been  used  to 
evaluate  the  technique.  The  apparent  fracture  toughness  of  this  material 
increases  with  increasing  crack  length  (R-curve  behavior).  This  behavior 
is  analysed  using  established  equations  in  conjunction  with  a  new  graphi¬ 
cal  method  that  obviates  the  need  for  measuring  the  apparent  toughness  of 
samples  containing  very  large  cracks.  The  fracture  toughness  at  "infinite" 
crack  length,  0.87  ±  0.08  MPa-mV2>  is  in  excellent  agreement  with  the  val¬ 
ues  obtained  by  conventional  methods.  The  effect  of  sample  size  on  the  ac¬ 
curacy  of  the  results  is  thoroughly  discussed.  The  thickness,  t,  of  the  sam¬ 
ple  should  exceed  by  a  factor  of  ten  the  initial  crack  length  produced  by  the 
indentation.  Also,  since  samples  are  tested  at  various  indentation  loads,  P, 
the  ratio  of  their  thicknesses  should  be  chosen  to  satisfy  the  condition  t  « 
p2/3  The  importance  of  several  parameters  in  the  design  of  the  apparatus 
is  also  discussed. 


I.  INTRODUCTION 


Many  techniques  have  been  developed  for  measuring  the  fracture 
toughness  of  ceramic  materials,  the  most  common  ones  being  the  double 
torsion  technique  and  the  double  cantilever  beam  test.  These  methods  re¬ 
quire  the  use  of  relatively  large  samples,  and  in  this  sense  are  material  in¬ 
tensive.  It  often  happens  m  materials  research  and  development  programs 
that  the  quantity  cf  material  available  is  limited,  hence  a  method  using 
small  specimens  to  measure  fracture  toughness  would  be  preferred.  The 
indentation  toughness  method*  has  evolved  as  a  particularly  useful  one  be¬ 
cause  it  is  fast,  simple  and  can  provide  useful  data  from  quite  small  sam¬ 
ples.  Nevertheless,  there  are  limitations  associated  with  this  method, 
which  arise  mainly  from  uncertainties  in  the  measurement  crack  length, 
the  development  of  the  pattern  of  cracks  and  the  presence  of  residual 
stresses  at  the  surface.  Also,  the  indentation  toughness  method  cannot 
provide  information  on  the  increase  in  the  resistance  to  crack  propagation 
with  increasing  crack  length  (R-curve  behavior). 

Over  the  past  ten  years  or  so  several  different  kinds  of  tests  have  been 
developed  specifically  to  investigate  the  mechanical  behavior  of  small  spec¬ 
imens.  The  miniaturized  disk-bend  test  (MDBT),  is  but  one  of  several  such 
tests,  the  history  and  variety  of  which  have  been  recently  reviewed  by  Lu¬ 
cas2.  We  have  developed  our  own  MDBT  apparatus,  as  well  as  procedures 
for  measuring  the  yield  stress  of  metals  and  alloys3,  and  decided  to  explore 
the  possibility  of  utilizing  it  additionally  for  measuring  the  fracture  tough¬ 
ness  of  brittle  solids.  This  decision  was  motivated  by  our  recent  attempts  to 
develop  stronger  and  tougher  ZnS-based  ceramics  for  infrared  window  ap¬ 
plications4'6.  We  typically  prepare  small  quantities  of  ceramic  alloys,  e.  g. 
alloys  in  the  ZnS-Ga2S3  system,  and  use  the  indentation  toughness  method 
to  compare  the  fracture  toughness  of  these  alloys  as  a  function  of  composi¬ 
tion  and  the  processing  conditions.  The  information  so  gathered  is  quite 
useful,  but  limited  for  the  reasons  already  mentioned. 

We  therefore  decided  to  explore  the  possibility  of  using  the  "controlled 
flaw"  method  to  measure  the  fracture  toughness,  T»,  of  several  ceramic 
materials.  The  basic  procedure  in  conventional  testing  is  to  first  indent  the 
specimen  and  then  measure  its  macture  strength  in  a  bend  test,  such  as 
four-point  bending  or  ring-on-ring  bending.  The  latt°r  employs  circular 
specimens  and  has  been  used  to  measure  the  fracture  toughness  of  flat 
disks  of  glass  and  alumina7  and  yttria8.  The  equations  of  fracture  mechan¬ 
ics  are  then  appl:;d  and  T«,  is  extracted  from  analysis  of  the  data. 

In  this  paper  we  describe  our  adaptation  of  the  controlled  flaw  tech¬ 
nique  using  samples  scaled  down  in  size  for  testing  in  the  MDBT  appara¬ 
tus3.  The  test  is  actually  a  miniaturized  ring-on-ring  disk-bend  fracture 
toughness  test  (MDBFTT),  with  which  T„  can  be  evaluated.  The  samples 
used  in  this  test  are  disks  3  mm  in  diameter  and  ranging  .n  thickness  from 
300  to  700  jrm.  Fracture  toughness  can  be  obtained  by  this  method  without 
measuring  the  crack  lengths,  thereby  eliminating  the  uncertainties  associ¬ 
ated  with  such  measurements. 
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There  are  many  advantages  to  the  MDBFTT.  First,  it  provides  an  ef¬ 
ficient  method  for  testing  many  specimens  prepared  from  a  limited  supply 
of  material,  which  as  already  noted  is  extremely  valuable  in  a  materials 
synthesis  program  in  which  small  quantities  of  material  are  produced. 
Second,  when  the  size  of  the  ceramic  specimens  is  so  small,  the  indented 
crack  is  likely  to  be  the  predominant  flaw,  so  that  more  reliable  statistical 
results  can  be  obtained.  Third,  sample  preparation  is  easy,  and  no  special 
equipment  is  needed  beyond  that  required  for  making  samples  for 
transmission  electron  microscopy.  Additionally,  the  ring-on-ring  geometry 
in  the  miniaturized  test  described  here  has  the  same  advantages  of  this 
geometry  for  larger  samples,  namely  constant  stress  on  the  tensile  surface, 
easy  alignment,  more  accurate  measurement  of  the  stress,  and  a  biaxial 
state  of  stress,  which  obviates  the  need  for  a  special  orientation  of  the  crack 
(cf.  four-point  bending,  in  which  the  crack  plane  must  be  perpendicular  to 
the  tensile  axis). 

II.  EXPERIMENTAL 

The  MDBFTT  apparatus  is  a  minor  variation  of  the  test  facility  de¬ 
signed  in  our  laboratory  for  bend-testing  more  ductile  materials3;  it  is  illus¬ 
trated  schematically  in  Fig.  1.  The  load  is  transmitted  by  a  punch  (0.981 
mm  in  diameter)  and  a  1  mm  diameter  ball  bearing  through  a  stainless 
steel  loading  ring  with  an  outer  and  inner  diameters  of  0.9  and  0.66  mm, 
respectively.  In  the  original  test  rig  the  load  is  transmitted  to  the  sample 
directly  by  the  ball  bearing,  but  in  this  adaptation  it  serves  to  assure  the 
alignment  of  the  loading  ring  along  the  vertical  axis.  The  entire  loading 
assembly  is  directed  through  the  guide  cylinder,  which  is  bored  to  a  diame¬ 
ter  of  1.005  mm.  The  inner  diameter  of  the  lower  die  is  2  mm,  hence  the 
samples  are  supported  over  an  annulus  of  0.5  mm.  Both  the  die  and  the 
cylinder  are  made  of  D2  tool  steel.  Shims  of  different  thickness  are  placed 
between  the  upper  and  lower  housings  to  accommodate  specimens  thicker 
than  the  300  pm  recess  in  the  lower  die. 

The  material  used  for  these  initial  tests  was  ZnS  prepared  by  chemi¬ 
cal  vapor  deposition  (CVD)  and  obtained  from  Morton  International.  ZnS 
was  chosen  because  it  is  commercially  available,  provides  a  reference  ma¬ 
terial  in  our  research  on  the  development  of  ZnS-based  infrared  transmit¬ 
ting  ceramics4*6,  and  has  well-characterized  physical  and  mechanical 
properties,  including  its  fracture  toughness  as  measured  by  conventional 
methods.  Selected  properties  of  this  material  are  listed  in  Table  1. 

The  material  was  first  cut  into  slices  varying  from  350  to  800  pm  in 
thickness,  depending  on  the  ultimate  thickness  desired  for  the  disk  speci¬ 
mens,  using  a  low-speed  diamond  saw  (South  Bay  Technology  model  650). 
All  the  slices  were  cut  along  a  plane  which  is  perpendicular  to  the  CVD 
growth  plane,  thus  the  normal  direction  of  the  disks  is  perpendicular  to  the 
CVD  growth  direction.  There  is  no  difference  in  the  hardness  parallel  and 
perpendicular  to  the  growth  direction.  Disks  3  mm  in  diameter  were  then 
trepanned  from  the  slices  using  an  ultrasonic  disk  cutter  (Gatan,  model 
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601)  with  SiC  powder  as  an  abrasive  slurry.  The  disks  were  then  metallo- 
graphically  polished  through  0.05  pm  alumina  powder.  The  specimens 
were  examined  under  the  optical  microscope  after  final  polishing  to  assure 
that  the  surfaces  were  scratch-free  and  that  no  other  pre-existing  flaws 
were  present.  The  thicknesses  of  the  specimens  were  chosen  according  to 
the  magnitude  of  the  indentation  load;  300  pm  thick  samples  were  used  for 
the  smallest  load  (0.98  N),  whereas  for  the  maximum  load  (19.6  N)  it  was 
necessary  to  use  specimens  700  pm  thick.  The  thicknesses  of  the  samples 
were  measured  to  within  ±  2  pm  using  a  micrometer. 

Vickers  indentations  were  made  in  the  center  of  the  disk  using  a 
Buehler  Micromet  microhardness  tester  for  loads  up  to  1  kg,  and  a  Tukon 
microhardness  tester  for  larger  loads.  The  same  loading  time  was  used 
over  the  entire  range  of  loads,  from  100  g  to  2  kg  (0.98  to  19.6  N).  The  centers 
of  the  disks  were  located  to  within  ±  10  pm  before  indentation  using  the  mi¬ 
crometer  stage  built  into  the  hardness  tester.  Unindented  disks  were  also 
tested  to  measure  their  flexural  strength.  At  least  five  specimens  were 
tested  for  each  indentation  load.  All  the  specimens  were  examined  in  the 
optical  microscope  after  testing  to  check  whether  or  not  the  Vickers  inden¬ 
tations  served  as  the  origins  of  the  failures.  Results  from  samples  that  did 
not  meet  this  criterion  were  discarded  (roughly  20  %  of  the  tests  run  at  the 
smallest  indentation  loads,  100  and  200  g  fell  into  this  category,  but  at  loads 
of  500  g  or  more  nearly  95  %  of  the  tests  were  acceptable).  In  most  of  the 
successful  tests  the  specimens  fractured  into  four  pieces;  a  representative 
photograph  is  shown  in  Fig.  2.  Specimens  indented  with  loads  of  500  g  of 
more  had  cracks  emanating  from  the  comers  of  the  indentations.  This  is 
probably  responsible  for  the  greater  proportion  of  successful  tests  at  the 
higher  loads. 

The  disks  were  placed  on  the  lower  die  with  their  indented  side  sub¬ 
jected  to  biaxial  tensile  loading.  A  small  amount  of  vacuum  grease  was 
placed  on  the  loading  ring  and  lower  die  to  lubricate  the  areas  of  contact 
with  the  specimen.  This  minimizes  the  lateral  forces  on  the  samples  dur¬ 
ing  loading.  The  bending  tests  were  performed  as  quickly  as  practically 
possible  after  indentation,  normally  less  than  30  minutes.  A  table  model 
Instron  testing  machine  (Model  TM-S-L)  with  crosshead  speed  of  8.4  pm/s 
was  used  for  all  the  bend  tests.  The  applied  force  was  measured  by  a  45.4  kg 
load  cell  (Interface  model  SM-100)  with  an  accuracy  of  1.8  g  (0.01766  N),  and 
the  displacement  was  measured  to  within  0.06  pm  using  a  linear  variable 
differential  transformer  (LVDT).  The  signals  from  the  load  cell  and  the 
LVDT  were  collected  by  an  IBM  PC  with  a  data  acquisition  rate  of  20  data 
points  per  second. 

m.  RESULTS 

A.  Measurement  of  the  Fracture  Stress,  Of 

A  representative  curve  of  applied  load,  F,  vs.  displacement,  w,  is 
shown  in  Fig.  3a.  The  sample  in  this  case  was  600  pm  thick,  indented  with 
a  load,  P,  of  4  .98  N,  but  the  features  seen  here  are  independent  of  these  pa- 
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rameters.  The  curve  is  typical  in  that  an  initial  nonlinear  region  at  the 
smallest  loads  and  displacements  is  followed  by  another  region  which  is 
linear  until  the  sample  fails  catastrophically.  Immediately  after  failure  the 
load  drops  so  rapidly  that  the  data  acquisition  system  cannot  record  it.  The 
maximum  load,  Fm,  (Fig.  3a)  was  subsequently  used  to  calculate  the  frac¬ 
ture  stress.  Of. 

A  few  samples  experienced  an  intermediate  drop  in  load,  followed  by 
an  increase  to  even  higher  loads  prior  to  catastrophic  failure.  An  example 
of  this  type  of  behavior  is  shown  in  Fig.  3b.  We  were  unable  to  uncover  any 
systematic  reasons  for  this  behavior;  for  example,  most  of  the  specimens 
that  behaved  this  way  also  fractured  into  four  pieces,  as  in  Fig.  2.  Apart 
from  this  aberration,  if  the  samples  exhibiting  this  behavior  were  normal 
in  all  other  respects  there  was  no  reason  to  exclude  their  values  of  Fm, 
hence  they  were  retained.  It  is  not  common  practice  to  determine,  or  at 
least  report,  load-displacement  curves  in  fracture-toughness  investigations 
of  ceramics;  ws  are  unaware  of  any  published  examples.  It  is  possible  that 
the  intermediate  load  drop  is  associated  with  the  temporary  arrest  of  one  of 
the  two  cracks,  but  there  is  no  obvious  way  of  verifying  this  suggestion. 

The  fracture  stress  was  calculated  from  the  formula  for  a  thin  circu¬ 
lar  plate  of  thickness  t  subjected  to  elastic  bending  under  the  simply-sup¬ 
ported  edge  boundary  condition9,  namely 
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where  v  is  Poisson's  ratio,  and  a  and  b  are  the  radius  of  the  lower  die  and 
the  loading  ring,  respectively.  Equation  (1)  is  valid,  strictly  speaking,  for  a 
loading  ring  with  an  infinitely  thin  wall  thickness,  which  is  not  realized  in 
practice.  Therefore,  a  formula  has  been  derived  for  the  more  realistic  case 
of  a  loading  ring  with  inner  and  outer  radii  6i  and  £>2.  respectively,  namely 
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It  turns  out  that  if  6  in  equation  (1)  is  taken  as  ( b\  +  62)/2,  equations  (1)  and 
(2)  yield  almost  identical  results.  The  difference  between  them  is  less  than 
0.1%  for  the  test  parameters  in  our  experiments,  therefore  equation  (1)  was 
used  for  all  the  calculations,  with  b  =  (b\  +  b2)!2. 

B.  Evaluation  of  the  Fracture  Toughness 

All  the  data  on  Of  as  a  function  of  P  are  presented  in  the  log-log  plot  in 
Fig.  4,  where  it  is  evident  that  Of  decreases  with  increasing  P.  This  type  of 
behavior  is  representative  of  ceramics  exhibiting  R-curve  behavior.  The  ba¬ 
sis  of  R-curve  behavior  in  controlled  flaw  tests  is  that  the  total  stress  inten¬ 
sity  factor,  K,  driving  crack  propagation  can  be  written  as7 
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where  c  is  the  crack  length,  cra  is  the  applied  stress,  y  is  a  dimensionless 
constant  related  to  the  crack  geometry  and  %  is  an  elastic  parameter10 
which  is  given  by 

/E\1V2 

X  =  0.004  (gj  ,  (4) 

where  E  is  Young's  modulus  and  H  is  the  Vickers  hardness.  The  first 
term  in  equation  (3)  is  due  to  the  applied  stress  and  the  second  is  due  to  the 
residual  stress  associated  with  the  indentation.  The  description  of  R-curve 
behavior  is  completed  by  equating  (3)  to  the  fracture  resistance,  T,  which  is 
a  function  of  the  crack  length,  i  s.  T  =  T(c).  A  crack  can  propagate  in  an 
unstable  manner  as  soon  as  it  reaches  a  critical  length,  c*,  defined  by  the 
conditions 


K(c*)  = 

=  Tic*) ; 

(5a) 

/dK\ 
l^dc  Jc  _  c* 

(5b) 

Equation  (5a)  defines  T(c*),  the  apparent  fracture  toughness  of  a  material 
containing  a  crack  of  length  c*.  Unstable  crack  growth  in  this  material  oc¬ 
curs  when  the  first  term  in  equation  (5b)  is  greater  than  the  second. 

The  principal  feature  that  distinguishes  the  various  theoretical 
treatments  of  R-curve  behavior  is  the  function  chosen  to  represent  T(c). 
Perhaps  the  most  common  function  chosen  represents  power-law  behavior, 
namely11-13 

T(c)  =  TjfeJ ,  (6) 

where  T0,  c0  and  n  are  empirical  material  constants  with  no  special  physi¬ 
cal  significance  except  that  T0  and  c0  represent  some  characteristic  tough¬ 
ness  and  crack  length,  respectively.  A  different  relationship  between  T(c) 
and  c  has  been  suggested  by  Cook  et  al.14,  namely 

T(c)  =  T„-§|.  (7) 

where  the  product  pQ  represents  a  "microstructural"  driving  force  and  Too 
represents  the  toughness  at  infinite  crack  length.  Other  functional  depen¬ 
dencies  for  T(c)  have  been  proposed15,  but  they  will  not  be  considered  fur¬ 
ther  here. 


It  is  not  difficult  to  show  using  equations  (5)  (see  e.  g.  the  paper  by 
Krause11)  that  the  function  in  equation  (6)  leads  to  a  power  law  relationship 
between  cr f  and  P  of  the  form 


Of  =  AP'q  , 

where  q  and  n  are  related  by  the  formula11 

1  — 3q 
n  “  2  +  2q  ' 


(8) 


(9) 


Whereas  it  is  certainly  possible  to  fit  a  straight  line  to  the  data  in  Fig.  4,  we 
do  not  believe  that  the  fit  is  good  enough  to  justify  the  description  of  R-curve 
behavior  by  equation  (6),  and  have  instead  chosen  to  use  equation  (7)  for  fur¬ 
ther  analysis. 

Cook  et  al.14  divide  the  fracture  behavior  into  two  regimes,  a 
"microstructure-controlled"  regime  in  which  Of  is  independent  of  P  but  the 
apparent  fracture  toughness  is  not,  and  an  "indentation-controlled”  regime 
in  which  <Jf  varies  as  P'P3  while  the  fracture  toughness  is  independent  of  P; 
a  critical  load,  P*,  separates  the  two  regimes.  They  observed  that  the 
strength  of  numerous  ceramics  approached  the  microstructure-controlled 
regime  as  the  indentation  size  began  to  approach  the  grain  size  of  the  mate¬ 
rial.  Using  equation  (7)  in  conjunction  with  equations  (5),  they  showed  that 
Of  is  related  to  by  the  equation 

3  /Too  ^4/3  1 

~  MfyVZ  l  4  J  (P  +  P*)U3  ’  (10) 

where  xP*  =  |J-Q.  Equation  (10)  predicts  the  existence  of  two  regimes  of  be¬ 
havior  in  a  plot  of  log  Of  vs.  log  P;  at  small  indentation  loads  of  should  ap¬ 
proach  a  constant  value,  whereas  at  large  loads  Of  should  vary  as  P'P3. 

We  believe  that  the  data  in  Fig.  5  are  more  consistent  with  the  predic¬ 
tions  of  equation  (10)  than  equation  (8),  with  Of  approaching  a  constant 
value  at  small  P  and  approaching  the  predicted  p-^3  dependence  at  large  P. 
In  previous  work14  the  values  of  T«  and  P*  have  been  determined  by  opti¬ 
mizing  the  fit  between  equation  (10)  and  the  experimental  results.  How¬ 
ever,  on  rewriting  equation  (10)  as 

3-<Ejts),<p+p*>-  <n> 

it  is  apparent  that  the  data  should  fall  on  a  straight  line  if  plotted  in  the 
form  l/Of  vs.  P,  and  that  Too  and  P*  are  readily  obtained  from  the  slope  and 
the  intercept  of  the  curve.  This  method  of  evaluating  the  data  overcomes 
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the  objection  raised  by  Krause11  about  the  need  for  large-scale  crack  behav¬ 
ior  in  the  evaluation  of  T^. 

The  data  in  Fig.  4  are  replotted  in  this  way  in  Fig.  5,  where  the  line 
drawn  is  a  regression  curve  obtained  by  least-squares  fitting  of  the  data.  To 
extract  the  best  values  of  T»  and  P*  from  the  slope  and  intercept  of  the  line 
in  Fig  5,  it  is  necessary  to  choose  the  correct  value  of  y.  Estimating  y  is  not 
easy,  because  this  parameter  is  related  to  the  stress  intensity  at  the  crack 
tip,  which  for  our  small  samples  is  very  possibly  strongly  influenced  by  the 
geometry  of  the  samples  as  well  as  the  crack.  In  the  following  analysis  we 
use  the  value  y  =  1.24  in  our  calculations,  as  recommended  by  Cook  et  al10. 
This  choice  is  justified  later. 

Using  the  values  of  H  measured  from  the  indentations  of  the  disk 
specimens  and  y  =  1.24  in  equation  (3),  the  values  of  TM  and  P*  obtained 
from  the  slope  and  intercept  of  the  fitted  line  are  =  0.87  ±  0.08  MPa-m172 
and  P*  =  3.07  N  (these  values  were  also  used  to  calculate  the  curve  in  Fig. 

4).  The  great  advantage  of  plotting  and  analyzing  the  data  in  this  way  is 
that  it  is  not  necessary  to  make  measurements  at  high  loads  (  P  >  P*)  in 
order  to  obtain  Too.  This  is  particularly  useful  for  the  MDBFTT,  since  scal¬ 
ing  up  the  size  of  the  samples  is  unnecessary. 

The  value  of  Too  =  0.87  ±  0.08  MPa-m1'2  is  dose  to  that  provided  by  the 
supplier  (1.0  MPa-m^2),  determined  by  the  Vickers  indentation  technique. 
It  is  also  in  good  agreement  with  values  obtained  by  more  conventional 
methods.  Too  =  0.67  MPa-m1'2  by  four-point  bending17  and  1.0  MPa-mL'2 
obtained  using  the  double  torsion  technique18. 

IV.  DISCUSSION 

Having  demonstrated  that  the  MDBFTT  yields  data  that  can  be  ana¬ 
lyzed  to  extract  a  value  of  Too  that  compares  favorably  with  those  obtained  by 
other  methods,  it  is  important  to  address  the  major  issue  of  the  effect  of 
spedmen  size  with  respect  to  the  flaws  introduced  by  indentation  and  the 
geometry  of  the  test  apparatus. 

There  is  a  minimum  acceptable  thickness  of  the  spedmen  in  the 
MDBFTT.  Anstis  et  al.1  have  shown  that  the  thickness  of  the  disks  should 
be  much  larger  than  the  diameter  of  the  initial  crack  to  eliminate  surface 
effects.  Their  criterion  is  t  >  10  Cq>  where  c0  is  the  equilibrium  crack  length 
after  the  indentation.  Thus  t  =  10  Cq  defines  the  lower  limit  of  the  thickness. 
In  prindple  there  is  no  upper  limit  to  the  thickness  of  a  disk  specimen  for 
the  ring-on-ring  test,  since  the  portion  of  the  disk  inside  the  loading  ring  is 
in  a  state  of  pure  bending.  However,  as  the  thickness  of  the  spedmen  is  in¬ 
creased  the  contact  stress  at  the  loading  ring  becomes  very  large,  because 
the  load  needed  to  produce  the  fracture  stress  is  large.  As  a  practical  mat¬ 
ter  we  arbitrarily  choose  as  a  criterion  a  maximum  thickness,  tm,  such 
that  the  maximum  contact  stress,  oc  =  Fm/2jibAb,  where  Ab  is  the  wall 
thickness  of  the  loading  ring,  does  not  exceed  twice  the  fracture  stress. 

,  Then,  setting  oc  =  2cf,  with  Of  given  by  equation  (1)  and  t  =  tm,  we  find  on 
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substituting  the  parameters  relevant  to  this  work  that  tm  =  830  |xm.  This  is 
somewhat  greater  than  the  thickest  sample  we  tested  (700  pm). 

The  geometry  of  the  sample  also  affects  the  values  of  fracture  tough¬ 
ness  through  the  parameter  y,  which  influences  the  stress  intensity  factor. 
To  estimate  we  used  the  empirical  equation  of  Newman  and  Raju16  to  cal¬ 
culate  the  stress  intensity  factors  of  surface  cracks  as  a  function  of  speci¬ 
men  geometry.  The  ratios  c*/t  and  c*/b,  as  well  as  <p,  where  <p  is  the  angle 
between  the  surface  of  the  specimen  and  the  crack  front  (measured  from 
the  mid-point  of  the  crack),  all  affect  the  stress  intensity  factor*.  In  o in¬ 
fests  the  effect  of  c*/b  is  small  because  b  »  c*.  If  we  assume  that  the  crack 
retains  its  semicircular  shape  during  propagation  and  assign  a  single 
value  of  the  stress  intensity  factor  along  the  crack  front,  then  the  equations 
of  Newman  and  Raju  can  be  used  to  calculate  the  stess  intensity  factor  as  a 
unique  function  of  c*/t. 

The  necessary  calculations  were  performed  by  computer  to  calculate 
the  value  of  relevant  to  the  parameters  characteristice  of  our  tests.  The 
required  estimates  of  c*  were  calculated  using  the  formula14 


c*  = 


X(P+P*) 


2/3 


(12) 


in  conjunction  with  the  experimentally  determined  values  of  P*  and  T^. 

For  the  range  of  c*/t  in  our  tests,  namely  0.11  to  0.33,  the  corresponding 
values  of  \|/  decrease  from  1.25  to  1.16.  Since  the  smaller  values  of  c*/t, 
hence  the  larger  values  of  ij/,  obtain  at  the  smaller  loads,  it  is  possible  that 
there  is  a  small  contribution  of  sample  geometry  to  the  observed  R-curve 
behavior.  Its  contribution  would  not  appear  to  be  significant,  however,  be¬ 
cause  if  we  use  the  average  value  y  =  1.205  to  calculate  TU,  the  result  is  T«  = 
0.85  MPa-m1''2,  which  differs  by  only  3%  from  the  value  obtained  using  \j/  = 


1.24. 


For  conventional  fracture-toughness  tests  on  plate-shaped  specimens 
the  thickness  and  crack  width  are  identical,  and  the  thickness  should  be 
much  larger  than  the  size  of  the  plastic  zone  at  the  crack  tip  to  assure  plane 
strain  conditions.  In  normal  bend  tests  on  pre-indented  large  ceramic 
samples  the  initial  and  critical  crack  sizes  are  ordinarily  so  much  smaller 
than  the  thicknesses  of  the  specimens  that  adjusting  the  size  of  the  speci¬ 
mens  to  accommodate  the  larger  flaws  produced  by  high  indentation  loads 
is  unnecessary.  For  example,  in  Krause's  experiments11,  the  ratio  of  crack 
size  to  specimen  thickness  varied  by  about  an  order  of  magnitude  (0.02  to 
0.20),  yet  decreased  by  only  6.7%  (from  1.20  to  1.12).  In  the  MDBFTT  the 
dependence  of  xi/  on  c*/t  is  more  sensitive  than  this,  with  decreasing  by 


* 

The  stress  intensity  factor  along  the  crack  front  is  not  constant,  and  is 
higher  at  regions  of  the  crack  tip  near  the  specimen  surface.  This  will 
affect  the  shape  of  the  crack  as  discussed  by  Krause11. 
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7.2%  as  c*/t  increases  by  a  factor  of  three.  This  suggests  that  the  thick¬ 
nesses  of  the  samples  in  the  MDBFTT  have  to  be  chosen  rather  more  care¬ 
fully  than  in  bending  tests  on  large  samples,  and  that  a  criterion  for  mak¬ 
ing  the  correct  choices  would  be  useful. 

We  suggest  as  an  appropriate  criterion  that  the  ratio  of  critical  crack 
length  to  the  sample  thickness  should  be  maintained  nearly  constant  to 
eliminate  specimen  size  effects  on  the  experimental  value  of  T«,.  Since  the 
critical  crack  length,  c* ,  is  not  known  beforehand,  our  suggestion  would  not 
appear  to  be  very  helpful  in  designing  experiments**.  However,  it  has 
been  reported  that  c*  is  proportional  to  Co  and  the  ratio  of  these  two  lengths 
is  independent  of  the  indentation  load14.  Since  Cq  is  a  function  of  P  through 
the  established  equations  of  indentation  toughness  measurements1,  we  can 
write 

p2/3  =  pct,  (13) 


where  p  is  a  proportionality  constant.  To  maintain  the  ratio  cVt  constant 
for  two  specimens  with  indentation  loads  Pi  and  P2,  the  relation  between 
the  thicknesses  of  the  two  samples  should  be 


ti 

t2 


4 


1  /Plf/3 

rw  • 


(14) 


Equation  (14)  provides  the  criterion  we  seek  because  it  specifies  the  relation¬ 
ship  between  t  and  P  required  to  maintain  cVt  constant  without  needing  to 
know  anything  about  the  fracture  behavior  of  the  material.  This  criterion, 
in  conjunction  with  t  >  10  Cq,  should  serve  as  valuable  guidelines  for  the 
preparation  of  specimens  for  the  MDBFTT. 

Finally,  we  address  the  issue  of  the  geometrical  parameters  used  in 
the  test  apparatus,  because  these  can  also  influence  the  results.  The  ratio 
of  the  diameters  of  the  lower  die  (2a)  and  the  loading  ring  (26)  is  the  most 
important  factor  here.  If  6/ a  is  too  small  the  contact  stress  exerted  on  the 
specimen  by  the  loading  ring,  and  vice  versa,  can  exceed  the  yield  stresses 
of  either  the  loading  ring  or  the  specimen.  On  the  other  hand,  from  practi¬ 
cal  considerations  the  inner  diameter  of  the  loading  ring  should  be  larger 
than  the  length  of  the  initial  crack  length,  i.  e.  b  »  Cq.  This  ensures  that 
the  area  over  which  the  stress  is  constant  is  as  large  as  possible.  However, 
the  value  of  b  la  does  have  an  upper  limit  because  <3{  becomes  zero  as  6  la  ap¬ 
proaches  unity,  as  can  be  seen  from  equation  (1).  It  does  not  seem  possible 
to  determine  an  optimum  value  of 6 fa,  but  as  a  working  hyopothesis  we 


**The  crack  lengths  c*  and  c1  are  not  necessarily  identical.  The  former  is 
calculated,  in  our  case  using  equation  (12),  whereas  the  latter  is  intended 
to  represent  the  actual  critical  crack  length  in  the  specimen. 
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suggest  that  if  6 la  satisfies  the  condition  0.3  <  6/a  <  0.7  the  accuracy  and 
sensitivity  of  the  test  is  assured;  6/a  =  0.45  in  our  apparatus. 

V.  CONCLUSIONS 

1.  It  has  been  demonstrated  that  the  fracture  toughness  of  a  ceramic 
material  can  be  measured  successfully  by  the  MDBFTT.  Experiments  on 
CVD -grown  ZnS  yield  the  fracture  toughness  TM  =  0.87  ±  0.08  MPa-m1^2, 
which  is  in  good  agreement  with  the  values  obtained  by  the  bend  testing  of 
large  specimens,  and  by  the  indentation  toughness  method. 

2.  The  data  on  ZnS  exhibit  R-curve  behavior.  The  equations  of  Cook 
et  al.14  accurately  describe  the  relationship  between  fracture  stress  and 
indentation  load.  A  graphical  procedure  has  been  developed  to  analyze  the 
data  using  the  equations  of  Cook  et  al.,  which  has  the  advantage  that  data 
at  large  indentation  loads  are  not  required  to  extract  the  correct  value  of  Too. 

3.  A  criterion  for  selecting  the  thicknesses  of  samples  in  the 
MDBFTT  is  proposed.  As  the  P  increases,  t  should  be  increased  as  P2^. 
This  assures  that  the  stress  intensity  at  the  front  of  a  crack  of  critical  size  is 
independent  of  the  thickness  of  the  specimens.  As  demonstrated  by  our 
results,  small  departures  from  this  criterion  do  not  seriously  affect  the 
values  of  Too. 
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Table  1.  Physical  Properties  of  the  CVD-Grown  ZnS  used  in  this 
investigation. 


Grain  Size 
(pm) 

Hardness 

(GPa) 

Flexural 

Strength 

(MPa) 

Young’s 

Modulus 

(GPa) 

Poisson’s 

Ratio 

Fracture 

Toughness 

(MN-m1/2) 

20-35 

1.5 

50-60 

74.5 

0.28 

1.0 

Fig.  1.  Schematic  diagram  of  the  apparatus  used  for  the  miniaturized  disk 
bend  fracture  toughness  tests. 


Fig.  2.  Photograph  of  the  pieces  of  a  700  pm  thick  sample  of  ZnS  after  frac¬ 
ture.  The  indentation  load  was  19.6  N. 


Fig.  3  Curves  of  applied  load,  F,  vs.  displacement,  w,  for  two  specimens 
indented  with  a  load  of  4.9  N.  The  curve  in  (a),  for  a  sample  570  urn 
thick,  was  typical  of  that  obtained  for  most  of  the  specimens.  The 
curve  in  tb),  for  a  sample  560  pm  thick,  exhibited  intermediate  un¬ 
loading  before  frac:ure. 


Pig.  4  A  log-log  plot  of  fracture  stress,  cf,  vs.  indentation  load,  P.  The 

shaded  rectangle  indicates  the  values  of  <?f  for  the  unindented  spec 
imens.  The  value  of  P*  and  the  solid  curve  were  derived  from  anal 
ysis  of  the  data,  as  described  in  the  text.  The  slope  of  the  dashed 
Line  is  -1/3. 


Fig.  5  The  data  in  Fig.  4  replotted  in  the  form  (l/cfp  vs.  P.  The  shaded 

rectangle  at  P  -  0  represents  the  data  on  unindented  samples;  these 
were  not  included  ,n  the  least-squares  fit  used  to  obtain  the  straight 
line  shown. 


Experimental  and  Theoretical  Studies  of  Second-phase 
Scattering  in  IR  Transmitting  ZnS-based  Window 

William  W.  Chen,  B.  Dunn,  J.  M.  Zhang 

Department  of  Materials  Science  and  Engineering 
University  of  California,  Los  Angeles,  CA  90024,  U.  S.  A. 


ABSTRACT 

A  model  based  upon  Mie  theory  was  developed  to  calculate  the  level  of  scattering  caused  by  the 
incorporation  of  a  second-phase.  Various  transmission  curves  were  calculated  in  order  to  study  the 
microsu  jc  ural  scattering  effects  of  pores  and  ZnG^Sj  particles.  The  model  was  able  to  reproduce  the 

spectral  tr  nsmission  curve  (2.5  -  10  pm)  of  a  hot-pressed  ZnS  sample  containing  a  porous  scatterer  only 
when  a  size  !:  'ribution  effect  was  included  in  the  model  bimodal  distribution  of  pores  was  found  to 
be  responsible  or  he  detcc.ed  scattering. 


1.  INTRODUCTION 

The  design  of  a  multiphase  ZnS-based  ceramic  window  is  influenced  by  a  compromise  between 
mechanical  and  optical  considerations.  This  process  -eeks  to  op'i^izc  the  strength  and  fracrure 
toughness  of  the  window  while  minimizing  the  scattering  due  to  the  second-phase.  From  an  optical 
viewpoint,  a  single-ph-'se  material  with  a  cubic  crystal  structure  is  best  suited  as  a  window  material 
because  the  anisotropic  and  second-phase  scattering  mechanisms  are  intrinsicr  Uy  absent.  However  the 
motivation  to  develop  a  multiphase  or  composite  ZnS-based  window  is  to  improve  its  mechanic ll 
properties  beyond  those  cr  u.c  single-phase  material.  The  approach  for  the  design  of  such  a  window  and 
some  mechanical  property  measurements  has  been  discussed  previously.1  There  is  usually  an  optimum 
value  of  parameters  such  as  the  volume  fraction,  size  and  size  distribution  that  results  in  the  best 
mechanical  properties.  Unfortunately,  the  features  which  characterize  the  dispersion  of  the  second-phase 
are  not  necessarily  the  most  suitable  ones  for  optical  transmission. 

Since  the  optical  scattering  is  very  sensitive  to  the  characteristics  of  the  dispersion  of  scatterer,  it 
is  important  to  be  able  to  understand  the  influences  of  each  second-phase  parameter  upon  the  optical 
transmission  of  the  window.  The  present  study  has  used  Mie  theory  as  the  basis  for  an  optical  scattering 
model.  The  model  enables  the  calculation  of  the  spectral  in-line  transmission  of  2.5  -  10  pm  radiation 
through  a  ZnS  matrix  while  varying  the  characteristics  of  the  c  spersion  before  the  window  is  fabricated. 
The  application  of  this  model  will  help  determine  the  microstructuia!  features  that  minimize  the 
scattering  loss  from  a  second-phase. 


2.  THEORETICAL  BACKGROUND 

Figure  1  is  a  schematic  of  the  various  losses  in  intensity  suffered  by  an  incident  parallel  beam 
( I J  when  it  traverses  a  matrix  containing  a  dispersion  of  phencal  particles.  The  ratio  of  the  detected 
intensity  (I)  over  the  incident  intensity  is  the  in-line  transmission  (T)  and  is  described  by  the  following 
equation: 

T  --  7  =  0  -  R)2  expi-yt)  (1) 

'o 

where  R,  y.  and  t  arc  the  reflectivity,  attenuation  or  extinction  coefficient  and  thickness  of  the  sample, 
respectively.  The  pre-exponential  term  accounts  for  the  loss  due  to  two  fresnel  reflections  (one  entering 


and  one  exiting  the  slab)  for  an  uncoatcd  sample.  The  reflectivity  is  determined  from  the  refractive 
index  by  use  of  the  following  equation: 


R  = 


(n-J)t 


(2) 


(n  +  1)2 

The  wavelength  dependence  of  the  refractive  index  must  be  known  in  order  to  calculate  the  spectral 
reflectivity.  For  polycrystalline,  cubic  ZnS,  the  refractive  index  is  approximately  2.25  over  the 
wavelength  range  2.5  -  10  pm.2  By  the  substitution  of  the  value  2.25  into  equation  2  and  by  assuming 
that  the  extinction  coefficient  is  zero,  the  transmission  from  Equation  (1)  is  calculated  to  be  about  12%. 


The  extinction  coefficient  is  related  to  the  number  of  particles  per  unit  volume  (N^  in  the  matrix 

by:3 

Y^C^  (3) 

where  is  the  scattering  cross-sectional  area  of  a  single  particle.  When  this  quantity  is  divided  by  the 
geometric  area  (tit2).  the  dimensionless  quantity  Q^,  the  scattering  efficiency  coefficient,  is  obtained. 
The  extinction  coefficient,  defined  in  terms  of  is  given  by:3 

3  V 

Y=47Qsc*  W 

where  V  is  the  volume  fraction  of  the  second-phase.  An  examination  of  Equation  (4)  indicates  that  the 
spectral  curve  must  be  computed  in  order  to  calculate  the  transmission. 


Rayleigh4  derived  an  expression  for  the  calculation  of  in  terms  of  a  size  parameter  (a)  and 
the  relative  refractive  index  (n)  for  the  scattering  of  electromagnetic  radiation  by  small,  dielectric 
spheres.  This  scattering  is  based  upon  the  condition  that  the  electric  field  within  the  sphere  is  uniform 
and  parallel  to  the  external  field  when  illuminated  by  polarized  radiation.  If  the  incident  field  oscillates 
harmonically,  the  particle  can  be  considered  an  oscillating  dipole  that  radiates  scattered  waves  in  all 
direction.  For  a  Rayleigh  scatterer  dispersed  in  air,  the  scattering  efficiency  factor  is 


8  „  ✓  n2  -  1  -.2 


_  o  .  /  ir  -  j  w 

Qscar  =  a  a  (  _2  ,  i) 


(5) 


3~  Vn2+F 

The  size  parameter  is  described  by  (2xr/X),  and  the  relative  refractive  index  is  (n^nj)  where  n,and  n2  are 
the  particle  and  matrix  refractive  indices,  respectively.  Equation  (5),  however,  is  only  valid  when  the 
upper  limit  of  the  particle  size  to  wavelength  ratio  (r/X)  is  less  than  or  equal  to  0.05.  For  2  pm  radiation, 
particles  of  0.1  pm  radius  or  smaller  are  considered  Rayleigh  scatterers,  and  in  that  region,  the 
transmission  curve  exhibits  a  characteristic  X-4  dependence  derived  from  the  quantity  a4  in  Equation  (5). 


The  sizes  of  the  second-phase  particles  in  ceramic  bodies,  however,  are  not  always  small  with 
respect  to  the  incident  radiation.  The  scattering  of  an  electromagnetic  wave  by  a  homogeneous  sphere  of 
arbitrary  size  is  derived  from  the  solution  of  Maxwell's  equations  with  appropriate  boundary  conditions. 
Mic  had  derived  the  solution  for  this  type  of  scattering,  and  the  scattering  efficiency  coefficient  is  given 
as  an  infinite  series  of  Legendre  polynomials  whose  coefficients  are  complicated  functions  of  Bessel 
functions.  The  Mie  formula  gives:5 

Qscat  =  ^  £2k+  I)(iak|2  +  (bk|2)  (6) 

k=l 

where  ak  and  b^  are  the  Mie  coefficients  and  are  functions  of  the  relative  refractive  index  and  size 
parameter. 


The  scattering  efficiency  factor  may  be  calculated  for  given  values  of  n.  a.  and  X  by  adding  the 
temis  from  Equation  (6)  in  series.  The  basic  problem  involved  in  Mie  theory  is  the  evaluation  of  the 


coefficients  and  b^.  Many  mathematical  methods6  are  available  to  determine  or  approximate  the  Mie 

coefficients,  but  most  evaluations  arc  now  done  numerically  with  the  aid  of  a  computer,7'8  the  summation 
converges  to  the  value  after  the  addition  of  a  certain  number  of  terms.  Tor  this  study,  the  spectral 
scattering  efficiency  coefficient  is  calculated  numerically.  Sample  calculations  of  Qv.>f  were  compared 
with  prior  studies  involving  systems  using  air  and  water  as  the  matrix  material.  A  comparison  of  the 
calculated  values  of  the  present  study  and  the  literature  values9  indicated  no  discrepancies. 


3.  EXPERIMENTAL 

High  density  ZnS  disks  were  fabricated  by  uniaxial  hot-pressing.  The  details  of  the  procedure  is 
described  elsewhere.10  The  sample  was  optically  polished  using  a  final  polishing  medium  of  0.05  pm 
alumina.  The  density  of  the  sample  was  determined  by  the  water  immersion  method,  and  the  infrared 
transmission  data  was  collected  with  use  of  a  Perkin-Elmer  FTTR  spectrophotometer. 


4.  CALCULATED  TRANSMISSION  CURVES  FOR  SINGLE  SIZE  SCATTERER 

The  parameters  which  affect  the  scattering  characteristic  of  the  two-phase  window  are  evident 
from  equation  (4).  Two  factors  to  consider  are  the  thickness  of  the  sample  and  the  volume  fraction  of  the 
scatterer.  They  are  wavelength  independent  which  means  that  increasing  or  decreasing  these  quantities 
only  shifts  the  transmission  level  higher  or  lower,  but  does  not  change  the  shape  of  the  transmission 
curve. 


The  scattering  efficiency  coefficient  is  wavelength  dependent  and  is  a  function  of  the  size 
parameter  and  the  relative  refractive  index.  The  amount  of  scattering  depends  upon  the  index  mismatch 
between  the  matrix  and  the  scatterer  and  also  upon  the  relative  size  of  the  scatterer  with  respect  to  the 
wavelength  of  the  incident  radiation. 

Transmission  curves  were  calculated  in  order  to  observe  the  influence  of  each  parameter  upon 
the  scattering.  In  each  of  the  figures,  a  series  of  curves  are  plotted  which  represent  the  transmission  when 
varying  the  value  of  one  of  the  parameters.  The  effect  of  volume  fraction,  size  of  the  scatterer  and  the 
type  of  second-phase  are  examined  over  spectral  range  2.5  -  10  pm. 

Figures  2  and  3  show  the  calculated  transmission  curves  of  a  ZnS  matrix  containing  porosity  as  a 
second-phase.  A  pore  (n(=  1)  was  chosen  as  the  scatterer  because  it  is  a  very  common  microstructural 
feature  in  ceramic  bodies.  Figure  2  is  a  plot  of  the  transmission  of  99,  99.5,  99.95  and  99.99%  dense  2 
mm  thick  ZnS  slabs  containing  0.3  pm  pores.  The  calculated  curves  indicate  that  for  window 
applications,  bulk  ZnS  pieces  must  be  fabricated  so  that  the  volume  fraction  of  voids  is  reduced  to  at  least 
lXlO"4.  This  illustrates  the  reason  why  it  is  often  difficult  to  fabricate  optical  quality  windows  by 
conventional  sintering  or  hot-pressing  methods  since  a  very  small  amount  of  residual  porosity  will 
adversely  degrade  the  transmission.  99.99%  dense  ZnS  window  will  scalier  about  20%  at  2.5  pm.  Even 
an  increase  of  an  extremely  small  amount,  4x10^  volume  fraction,  of  the  porous  phase  to  99.95% 
density  will  significantly  decrease  transmission  in  the  2.5  -  5  pm  range  (15%  transmission  @  2.5  pm). 
To  determine  the  transmission  curves  as  function  of  pore  size,  the  density  and  sample  thickness  are  fixed 
at  99.5%  and  2  mm,  respectively.  Figure  3  shows  that  scattering  losses  are  much  lower  for  smaller  size 
pores  and  that  only  pores  smaller  than  1000  A  are  essentially  homogeneous  to  2.5  to  10  pm  radiation. 

A  porous  phase  in  a  ceramic  piece  is  an  inherently  undesirable  characteristic  because  it  acts  as  a 
stress  concentrator  thus  degrading  the  intrinsic  strength  of  the  sample.  On  the  other  hand,  the  controlled 
incorporation  of  ZnGa2S4  precipitated  into  a  fully  dense  ZnS  matrix  is  of  interest  due  to  the  possible 


improvement  of  the  mechanical  properties  of  the  two-phase  material  over  the  pure  one.10  The  modelling 
of  the  transmission  curve  for  hypothetical  ZnS-based  windows  containing  different  characteristic 
enrembles  of  ZnGa,S4  precipitates  would  facilitate  the  determination  of  the  critical  volume  fraction  and 
size  necessary  to  achieve  a  useful  transmission  property. 

There  is  no  refractive  index  data  for  ZnGa^S.^  in  the  literature,  however,  the  refractive  index  of 
an  analogous  compound,  a  ZiilnjS,^  single  crystal,  has  been  measured  from  0.6  to  1.25  pm.11  The 
dispersion  of  its  refractive  index  as  the  IR  region  is  approached  is  small,  and  this  study  has  extrapolated 
the  ZnIn2S4  dispersion  curve  in  the  2.5  -  10  pm  range.  An  examination  of  the  extrapolated  curve 
indicates  that  a  refractive  index  of  2.5  is  a  good  value  for  use  in  the  model.  This  value  represents  the 
worst  case  refractive  index  because  the  index  of  ZnGa,S4  should  be  smaller  since  the  polarizability  of  the 

Ga3+  ion  is  smaller  than  the  In3+  ion. 

Figure  4  is  a  plot  of  transmission  versus  wavelength  of  a  2  mm  thick  ZnS-based  ceramic 
containing  500  A  size  ZnGa,S4  precipitates  as  a  function  of  volume  fraction.  This  size  of  second  phase 
particles  is  not  uncommon  for  other  systems  which  are  exploited  for  their  precipitation  strengthening  and 
hardening  effects.12  The  scattering  increases  as  the  volume  fraction  of  the  second-phase  increases, 
however,  most  of  the  scattering  is  limited  to  the  2.5  to  4  pm  range,  and  even  with  5%  volume  fraction  of 
precipitates  the  transmission  is  only  reduced  to  64%  at  2.5  pm.  The  transmission  curves  for  a  2  mm  thick 
ZnS  ceramic  matrix  containing  fixed  5%  volume  fraction  of  different  size  zinc  thiogallate  are  plotted 
versus  wavelength  in  figure  5.  Similar  manner  to  the  behavior  of  pore  scatterer,  the  transmission 
decreases  as  the  size  of  the  ZnGa2S4  precipitate  increases.  The  plot  reveals  that  5%  volume  fraction  of 
250  A  radius  particles  is  virtually  transparent  to  2.5  to  10  pm  radiation.  A  comparison  of  the 
hypothetical  0.1pm  size  zinc  thiogallate  and  pore  (Figures  3  and  5)  indicates  that  the  porous  sample 
scatters  less  despite  the  fact  that  the  absolute  value  of  the  refractive  index  mismatch  (An)  for  the  pore  is 
much  greater  than  that  of  the  ZnGa^.  We  also  see  that  the  number  of  scatterers  is  the  dominant 
influence  upon  the  scattering  characteristic  of  these  two  transmission  curves  because  the  porous  sample 
contains  an  order  of  magnitude  fewer  scatterers. 

There  are  still  many  combinations  of  relative  refractive  index,  volume  fraction  and  size  of  the 
scatterer  for  which  the  transmission  curves  have  yet  to  be  calculated.  By  the  use  of  the  model,  we  have 
seen  that  with  all  else  being  equal,  the  scattering  decreases  as  the  refractive  index  mismatch,  the  volume 
fraction  second-phase,  and  the  size  of  the  scatterer  decreases.  However,  with  the  proper  combination  of 
each  parameter,  it  is  possible  to  tailor  a  two-phase  microstructure  for  ZnS  which  is  transparent  to  2.5  to 
10  pm  radiation. 


5.  EXPERIMENTAL  RESULTS  AND  ANALYSIS 
5.1.  Fully  dense  samples  and  single  pore  size  scattering 

The  infrared  transmission  for  a  2.1  mm  thick  hot-pressed  ZnS  and  of  IRTRAN  2,  a 
commercially  available  hot-pressed  ZnS  window,  are  plotted  in  figure  6.  The  IRTRAN  2  sample 
exhibits  near  theo  etical  transmission  accounting  for  two  fresnel  reflections  at  the  surfaces.  This 
behavior  indicates  that  y  is  zero  and  that  the  material  is  fully  dense  (>  99.99%).  The  measured  densities 
of  the  hot-pressed  ZnS  (30.000  psi  pressure)  and  the  IRTRAN  2  samples  are  99.74  and  99.77%, 
respectively,  based  upon  the  x-ray  density  whose  lattice  parameter  is  5.406  A.13  The  corrected  density 
of  our  hot-pressed  sample  is  99.964%.  Addition  of  the  correction  factor  to  the  IRTRAN  2  measured 
density  would  make  its  density  consistent  with  the  infrared  transmission  data. 


The  transmission  data  for  the  hot-pressed  ZnS  sample  (30,000  psi  pressure)  was  analyzed  by  the 
scattering  model  in  an  attempted  to  derive  the  size  of  the  pore  and  thus  test  the  validity  of  the  model. 
However,  it  was  not  possible,  based  upon  the  model  and  the  experimental  data  to  determine  a  pore  size 
since  none  of  the  calculated  transmission  curves  fitted  well  with  the  measured  one.  This  inconsistency 
revealed  that  the  assumption  of  a  single  size  scatterer  for  real  systems  was  invalid,  and  therefore  a 
function  describing  the  size  distribution  of  the  scatterer  needs  to  be  incorporated  into  the  model. 

S.2.  Polydispersity  modification  and  experimental  curve  fitting 

The  number  of  particles  per  unit  volume  is  not  a  constant  anymore  but  is  in  now  a  function  of 
the  radius  (N(r)  =  NQ  f(r)).  Modification  of  Equation  (4)  based  upon  this  new  assumption  gives: 

y.M/w  Qsca  ^dr/Aw^dr}  (7) 

ro  ro 

where  r,  and  rQ  are  the  upper  and  lower  limits  of  the  size  distribution  function  (f(r)).  The  extinction 
coefficient  now  describes  the  scattering  due  to  an  ensemble  of  scattercrs  whose  size  is  described  by  f(r). 

An  iterative  process  was  used  to  fit  the  experimental  transmission  data  to  the  one  generated  by 
the  model.  The  f(r)  used  for  the  fitting  was  a  zeroth-order  lognormal  distribution  (ZOLD)  function  that  is 
characterized  by  two  parameters,  the  modal  radius  and  width  (r  and  a).13  It  was  found  that  a  bimodal 
ZOLD  distribution  generated  a  transmission  curve  illustrated  in  figure  7  using  the  appropriate 
experimental  sample  thickness  and  volume  fraction.  The  bimcdal  distribution  is  plotted  semi- 
logarithmically  in  figure  8.  There  is  rather  good  agreement  between  the  calculated  curve  and  the 
experimental  transmission  data  points. 


6.  SUMMARY 

An  optical  model  based  upon  Mie  theory  was  developed  to  study  the  microstructural  influences 
upon  the  2.5  -  10  pm  infrared  transmission  of  ZnS-based  ceramics.  It  was  found  that  the  scattering  losses 
are  minimized  as  the  refractive  index  mismatch,  volume  fraction  second-phase  and  size  of  uae  scatterer  is 
decreased,  however,  the  model  based  upon  the  use  of  a  single  size  scatterer  was  unable  to  calculate  a 
transmission  curve  which  fitted  with  the  experimental  transmission  curve  of  a  ZnS  sample  containing 
porous  scatterer.  A  model  which  incorporates  the  size  dispersion  effect  of  the  scatterer  was  able  to 
generate  a  transmission  curve  that  matches  the  experimental  results.  The  agreement  between  the 
experimental  and  calculated  curves  indicates  that  it  is  now  possible  model  to  predict  the  transmission 
behavior  for  ZnS-based  ceramics  prior  to  fabrication  of  the  window. 
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Fig.  1 .  Schematic  of  loss  mechanisms  upon  an  incident 
radiation  of  intensity  Iq  where  R  and  I  are  the  fresnel 
reflection  loss  and  the  transmitted  in-line  intensity, 
respectively 
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Fig.  2.  Transmission  curves  of  hypothetical  2  mm  thick  ZnS 
parallel  slab  containing  0.3  (pm)  pore  as  a  function  of  volume 
fraction  porosity 
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Fig.  5.  Transmission  curves  of  a  hypothetical  2  mm  thick 
ZnS  parallel  slab  containing  a  fix  5%  volume  fraction  of 
different  size  ZnGa2S4  particles 
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Fig.  6.  Experimental  transmission  curves  of  Irtran  2  and  of 
a  2.1  mm  thick  hot-pressed  ZnS  disk 


transmission 


Fig.  7.  The  experimental  transmission  curve  of  2.1  mm  thick 
ZnS  disk  hot-pressed  at  30,000  psi  is  indicated  by  the  square 
data  points.  The  solid  line  represent  a  transmission  curve 
generated  to  fit  the  experimental  data  based  upon  a  pore  distribution 
plotted  in  fig.  8. 


Fig.  8.  Semi-logarithmic  plot  of  a  bimodal  ZOLD  function 
versus  pore  radius.  The  modal  radii  are  0.08  and  1 .8  pm,  and 
the  widths  of  the  distribution  (a)  are  0.38  and  0.4  respectively. 


